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TREFACE 

I 

■ This  interim  report  describes  work  performed  by  personnel  of  the  University 

_ of  Missouri  Research  Reactor  Facility,  Columbia,  Missouri  65201  during  the 

period  from  30  June  1976  to  31  December  1976,  under  Contract  F33615-76-C-5230, 
Project  7371.  The  program  was  monitored  by  Dr.  Robert  J.  Spry,  Air  Force 
Materials  Laboratory,  Air  Force  Systems  Command,  Wright-Patterson  Air  Force 
Base,  Ohio. 

The  program  was  directed  toward  the  compensation  of  residual  boron  in 
infrared  detector  grade  silicon  by  nuclear  transmutation  doping.  This  in- 
vestigation was  conducted  by  Dr.  J.  M.  Meese,  principal,  investigator.  Dr. 

D.  M.  Alger,  Mr.  S.  L.  Gunn,  Mr.  W.  F.  Richardson,  Mr.  P.  J.  Glairon,  Dr. 

D.  M.  McKown,  and  Dr.  J.  Steven  Morris. 
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I.  ^INTRODUCTION 

The  primary  objective  of  this  program  is  to  devise,  through  the  research 
to  be  discussed,  an  optimum  set  of  techniques  capable  of  compensating  the 
residual  boron  always  found  in  detector  grade  silicon  by  the  process  of  nuclear 
transmutation  doping.  The  material  to  be  developed  by  this  program  is  intended 
for  use  in  infrared  detector  systems  such  as  those  of  the  Air  Force  LADIR  (Low 
Cost  Arrays  for  the  Detection  of  Infrared)  program. 

The  best  silicon  available  today  contains  residual  boron,  a shallow 
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acceptor,  in  concentrations  of  the  order  of  10  B/cm  . This  is  because  the 
segregation  coefficient  of  this  impurity  is  close  to  unity  and,  therefore,  will 
not  be  substantually  removed  by  successive  passes  in  a zone  refiner.  This 
residual  boron  concentration  has  several  adverse  effects  on  the  operation  of 
extrinsic  photoconductor  IR  detectors  fabricated  from  Ga  or  In  doped  silicon. 
Since  the  B ionization  energy  is  less  than  those  of  Ga  or  In,  the  B levels 
thermally  ionize  at  temperatures  below  those  impurities  of  interest  forcing 
the  detectivity  out  of  the  background  limited  condition  at  prematurely  low 
temperatures.  This  effect  requires  these  detectors  to  be  cooled  by  liquid 
helium  rather  than  by  a helium  gas  closed  cycle  refrigerator.  As  a result, 
system  size  and  cost  increase  unless  these  boron  levels  can  be  populated  by 
carriers  from  compensating  donor  levels.  Furthermore,  the  shallow  boron 
optical  transition  is  beyond  the  atmospheric  transmission  windows  through  which 
the  IR  data  information  is  obtained.  Therefore,  these  boron  optical  transitions 
can  respond  to  undesirable  background  thermal  radiation  to  produce  a type  of 
nonuseful  background  noise. 

Both  of  these  problems  can  be  eliminated  with  a minimum  degradation  in 
theoretical  detector  gain  if  nearly  exact  compensation  of  the  boron  acceptors 
is  accomplished  by  the  addition  of  shallow  donors.  Greatly  overcompensating 
the  boron,  however,  is  undersirable  because  the  additional  donor  carriers  fill 
the  Ga  or  In  levels  once  all  the  B levels  are  compensated.  This  reduces 
detector  gain  by  reducing  the  number  of  desirable  optical  transitions  from  which 
signal  information  is  gained.  Additionally,  the  extra  impurities  from  over- 
compensation represent  additional  charged  impurity  scattering  centers  reducing 
the  carrier  mobility  which  is  dominated  by  charged  impurity  scattering  at  the 
low  temperatures  at  which  these  detectors  operate.  This  lowering  of  carrier 
mobility  reduces  the  gain  of  the  detector  which  is  proportional  to  the  product 
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of  mobility  times  lifetime. 

Nuclear  transmutation  doping  of  silicon  appears  to  offer  several  advantages 
over  conventional  doping  techniques  for  producing  the  desired  compensating 
donor  concentration.  Briefly,  the  transmutation  doping  process  is  based  on  the 
conversion  of  the  isotope  ^Si  to  ^Si  through  an  (n,  y)  reaction  induced  by 


thermal  neutron  capture  during  reactor  irradiation.  The 
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stable  and  decays  through  6 emission  to 

31 


Si  isotope  is  un- 

P,  the  desired  donor  impurity.  Since 

"3T  1 

the  production  rate  of 
slow  enough  to  provide  careful  control  of  the  dopant  concentration.  Furthermore, 
the  donors  added  by  transmutation  doping  are  distributed  uni formally  through- 
out the  lattice  because  the  microscopic  abundance  ratios  of  the  Si  isotopes 
are  uniform  and  the  thermal  neutron  diffusion  length  is  very  long  compared  to 
typical  ingot  dimensions. 

The  two  major  advantages  of  NTD-Si,  tight  control  over  doping  concentration 
and  uniformity  of  dopant  distribution,  are  somewhat  offset  by  the  radiation 
damage  effects  which  occur  during  irradiation.  These  effects  cause  no  difficulty 
for  float  zone  Si  neutron  doped  to  a resistivity  of  100  ft-cm  or  lower.  For 
higher  final  resistivity  material,  the  success  or  failure  depends  intimately 
on  residual  impurity  defect  interactions  and  their  annealing  properties.  One  of 
the  major  goals  of  the  present  research  is  to  define  the  limits  over  final 
resistivity  which  can  be  achieved  relative  to  the  properties  of  the  starting 
material. 

The  balance  of  this  report  is  contained  in  five  major  sections  as  follows: 

• Section  II  develops  the  theoretical  tools  to  understand  the  expected 
electrical  properties  resulting  from  transmutation  doped  Si  as  a function 
of  neutron  fluence.  A discussion  of  the  transmutation  doping  process 

is  also  presented. 

• Section  III  is  a description  of  the  experimental  techniques  and  apparatus 
employed  throughout  this  work. 

• Section  IV  gives  a discription  of  the  experimental  results  and  a discussion 
of  possible  interpretations. 

• Section  V discusses  the  significant  accomplishments  achieved  during  the 
program  and  discusses  areas  of  difficulty  where  further  research  will  be 
required. 


II. 


THEORETICAL  ANALYSIS 


The  basic  physics  of  transmutation  doping  in  Si  has  been  discussed  by 
Tanenbaum  and  Mills.1  Although  they  discuss  all  possible  nuclear  reactions 
with  all  the  silicon  and  oxygen  isotopes  such  as  (n,  y),  (n,  2n),  (n,  p),  and 
(n,  a),  we  will  confine  our  attention  to  the  (n,  y)  reactions  with  Si  as  being 
the  only  significant  reactions  in  a thermal  neutron  flux.  This  viewpoint  has 
been  verified  experimentally.1 

The  natural  abundance,  thermal  neutron  absorption  cross  section  (n,  y) 
and  half-lifes  of  the  relevant  isotopes  of  silicon  and  phosphorus  are  listed 
in  Table  1. 

TABLE  1.  NUCLEAR  DATA 


Natural 

Abundance 

Isotope 

cr(n,  y) 

Tl/2 

1 

92.21% 

^Si 

0.08  + 0.03  b 

j 

stable 

k.68% 

29Si 

0.28  + 0.09  b 

stable 

3.05% 

3°Si 

0.11  + 0.01  b 

2.62h 

100% 

31P 

0.20  + 0.02  b 

lU.3d 

At.  Wgt . Si  = 28.086,  N(Si) 

-5 

Density  Si  = 2.33  gm/cm 
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5 x 10  atoms /cm 


The  nuclear  reactions  of  interest  are  the  following: 

28Si  (n,  y)  29Si  (1) 

29Si  (n,  y)  30Si  (2) 

30Si  (n,  y)  31Si *31P  + 6“  (3) 

T1/2  = 2.62h 

31P  (n,  y)  32P  ►32S  + 0“  (U) 

Tl/2  = lU*3d 

The  first  two  reactions  redistribute  the  natural  abundance  of  the  silicon 
isotopes  but  do  not  lead  to  any  doping  effects  or  radioactivity.  Reaction  (2) 


tmsmmst  igr; 
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is  quite  fortunate  in  that  it  serves  to  replace  the  Si  isotope  approximately 

31 

as  it  is  burned  up  in  reaction  (3)  which  produces  the  dopant  P.  Reaction  (4) 

31  32  - 

transmutes  the  stable  isotope  P to  the  unstable  isotope  P which  8 decays 

32 

to  S.  Since  reaction  (4)  is  a second  order  effect,  the  doping  is  thought 

31 

to  be  negligible  compared  to  the  doping  from  P (see  Ref.  l).  This  reaction 

is  important,  however,  as  the  primary  source  of  radioactivity  for  times  which 

are  long  compared  to  the  half-life  of  reaction  (3). 

From  the  data  given  in  Table  1,  we  can  calculate  the  production  rate  of 

3 

the  desired  dopant.  The  total  number  of  Si  atoms/cm  is  found  to  be 


op 

N(Si)  = 2.33  (6.023  x 10' ) 

28.086 

30, 
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4.996  x 10  atoms/cm  . 


The  density  of  Si  isotopes  is  then  found  to  be 

N(3°Si)  = 0.0305  N(Si)  = 1.524  x 1021  atoms/cm3. 

The  theoretical  production  rate  for  phosphorus  is  then  found  from  the 
expression 

d[31p]  = N( 3°Si ) o[30Si  (n,  y)  31Si] 

d$ 

= 1.676  x 10  ^ 3^P  atoms /cm3/n/cm2 

2 

where  $ is  the  irradiation  fluence  of  thermal  neutrons  (n/cm  ).  Table  2 lists 
the  theoretical  production  of  donors  in  Si  in  several  useful  sets  of  units. 


TABLE  2. 
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P THEORETICAL  PRODUCTION  RATES 


Production  Rate  (K) 

Units 

1.676  x 10~U 

1.676  x 1011+ 
3.355 

3 2 

P atoms /cm  per  n/cm 

P at-oms/cm3  per  lO^n/cm2 

ppb  per  10  n/cm 

Knowing  the  concentration  of  donors  added  per  unit  fluence,  it  is  now 
possible  to  calculate  the  final  resistivity  after  irradiation  and  annealing 
from  the  starting  resistivity  if  the  mobility  of  the  majority  carrier  is  known. 
To  this  end  an  analysis  of  mobility  data  for  Si  has  been  performed  on  data  taken 
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O o 

from  two  review  papers.  * Since  these  authors  quote  the  same  references  but 

give  different  values,  we  have  averaged  when  this  has  occurred.  We  have  also 

discarded  three  values  which  differed  abnormally  from  the  mean,  y^  = 13^5, 

y_  = 1590  and  vr  = 356  cm^/V-sec. 

Dn  Dp 

The  result  of  the  average  of  all  other  data  is  summarized  in  Table  3. 
TABLE  3.  TYPICAL  EXPERIMENTAL  MOBILITY  VALUES  IN  Si 


= ( 49h  +_  ll)  cm2'/V-sec 
y__  = (1396  + 6I4)  cm2/V-sec 

rp  = = 0.77  + 0.06 

r = 1.3  + 0.1 

2.826 

= (1.265  +_  0.29)  x 10  cm-0 

P 

= (253  + 5.7)  ppb 

P 

= (U.U77  + 0.197)  x 1015  cm-3 

P 

= (89.5^  ± 3.94)  ppb 

P 

These  expressions  will  be  useful  expressions  for  estimating  the  ppb  added  by 
neutron  irradiation  in  the  sections  to  follow. 

If  the  initial  and  final  carrier  concentrations  are  large  at 
room  temperature  compared  to  the  number  of  intrinsic  carriers,  then  the  following 
expressions  for  the  carrier  densities  in  terms  of  donor  and  acceptor  concentrations 
are 


y /y 
n p 


-1 


[ey  ] 
= L ^n 


-1 
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p-type  starting  material 


then. 


Po  = NA  - ND 

o 

n = N_  - Na 
o D H 


pf  = na  - ND  * 


nf  = NDf  " V 


nf  " no  = ND,  “ ND  * 
f o 


nf  + p0  = V - ND  * 
f o 


p0  - Pf  = N - Nd  , 
i o 


n-type  starting  material 
p-type  ending  material 
n-type  ending  material. 


n-type  to  n-type 
p-type  to  n-type 
p-type  to  p-type. 


But  Np  - ND  is  just  the  concentration  of  donors  added  by  transmutation  which 
we  P ° will  call  N_  = K $ where  K is  the  production  rate  given  in 
Table  2 and  <5  is  the  fluence  (number  of  neutrons  /cm  ).  Therefore,  we  see  that 
in  all  cases  for  which  n^  can  be  neglected,  the  concentration  of  phosphorus 
added  can  be  determined  by  the  concentrations  of  uncompensated  carriers  and  the 
material  type. 

We  can  determine  the  concentration,  Np,  of  phosphorus  atoms  added  either 
from  the  fluence  4>  or  the  initial  and  final  resistivities  and  types  using  the 
following  expressions, 

Np  = K $ = 3.355  x 10“18  SSk.  g $ (5) 

n/ cm 


K 

n 


K 

n 


n-type 


n-type 


p-type 


n-type 


(6) 


p-type 


p-type 
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where  K =89.5^  ppb  fi-cm  and  K = 253  ppb  fi-cm.  The  conversion  between  ppb 
n o ,P  3 

and  atoms/cin  is  1 ppb  = 5 x 101  atoms/cm  . 

By  equating  the  two  expressions  above,  the  fluence  necessary  to  obtain 

a desired  final  resistivity  can  be  determined  from  the  initial  resistivity  and 


type  as 


K 


K -1 

n 

K nf 


If  <Pp 


n-type 


p-type 


p-type 


n-type 


n-type 


p-type 


The  resistivity  after  irradiation  and  annealing  as  a function  of  fluence 
can  also  be  calculated  for  those  cases  where  n^  is  not  negligible  compared  to 
the  majority  carrier  concentration. 

The  resistivity  vs  fluence  has  been  calculated  for  transmutation  doped 
samples  using  the  following  assumptions: 

1.  100#  electrical  activity  is  obtained  by  annealing. 

2.  The  mobilities  remain  constant. 

2 2 

(hn  = 1500  cm  /V-sec,  = 600  cm  /V-sec)  in  the  following  figures. 

3.  The  shallow  donors  and  acceptors  are  completely  ionized  at  room  temperature. 

i 8 2 

k.  The  production  rate  of  phosphorus  is  given  by  K = 3.355  ppb/10  n/cm  , 

which  is  obtained  from  the  cross-section. 

5.  The  number  of  intrinsic  carriers  at  room  temperature  is 

10  -3 

ni  = 1.391  x 10  cm  . 

From  charge  conservation,  the  complete  ionization  assumption,  and  the 


relation 


we  obtain 


np  = nt 


p = 1/2[NA  - Nd  +v/(Na-Nd)2  + Hn.2] 


For  a production  rate  of  donors  given  by  K,  this  equation  becomes  for  p-type 


material 


p = 1/2[(Na  - Nd  - K4>)  + V(Na  - Nd  - K<f>)2  + 1»  n.2]  . 
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If  we  define  the  degree  of  initial  compensation  kQ  by  the  expression 

, N_ 
k = D , 

° na 

then  the  degree  of  compensation  as  a function  of  fluence  is  given  by 
k = ND  +K<f)  . 


NA 

For  complete  compensation,  k = 1;  for  no  compensation  k = o.  We  define  one 
more  quantity,  a,  the  degree  of  compensation  remaining  by 


a = !-k  = V ND  - ** 


N. 


N . 


For  no  irradiation,  we  assume  that  p >>  n.  so  that 

a - 1-*  - »A  - -D  0 - »i 

o o 


N, 


N, 


"A  A 

It  can  be  shown  by  substitution  in  Eq.  (l)  that 


Na 

P = A 

2 


a + 


/ — 7 — 

1 2 J 2n.  X 
va  ■*/__ i_ | 

VW  J 


(9) 


(10) 


We  now  define  the  fluence  necessary  for  complete  compensation  as  <f>  . This 
is  obtained  when  k = 1 or  a = o and  is  obtained  from  Eq.  (9) 

o = l-  l=  N -N  -K4> 

A D c 


or 


4>  = Na  - _ p 

rc  A D = *o 


It  can  now  be  shown  that 
N 


= (1  - WD)  (1  - <P  ) 


(i  - i ) • 


Then,  from  Eq.  (10) 


P = Vo  fl-iV  \/l  - <t»  )2  1 

2 IV  <J>C  \NAa|  J 


But  from  Eq.  (9)  with  <)>  = o, 

NAa  = p 
A o o 


i 


t 


! 


■I"  ■■  - ' * [ * 


(11) 


therefore,  we  obtain 


It  can  similarly  be  shown  for  overcompensation  (n-type  produced) 

n^r-  U-i)+  vXl  - (t>)2  + /2nl  \ 

2 L *C  *c  vt'l  -I 


(12) 


It  should  be  noted  that  the  number  of  carriers  is  symmetrical  with  respect 

to  the  fluence  on  either  side  of  d>  = d>  . It  should  also  be  noted  that  the 

c 

number  of  free  carriers  as  a function  of  fluence  is  independent  of  the  degree 

of  compensation.  It  can  be  shown  that  Eqs.  (ll)  and  (12)  form  a universal  doping 

curve  since  doubling  the  initial  concentration  doubles  the  normalized  flux 

coordinate  (l  - ),  i.e.,  if  p = 2 p then 

c o o 


P 


An  expression  for  the  resistivity  as  a function  of  fluence  can  be  obtained 
from  the  expression 

P = 1 

e(nyn  + pyp) 


and 

therefore , 


P = 


b « un/v 


P 


L 

/,  2 2. 

ey  (bn  + p ) 


ey  (bn2  + n.2) 
P i 


p-type 


n-type 


(13) 
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By  substituting  n or  p as  a function  of  fluence  from  Eqs.  (ll)  and  (12)  the 
resistivity,  p,  can  be  obtained  as  a function  of  fluence.  This  expression  is 
plotted  in  Figures  (l)  and  (2)  for  typical  p-type  float  zone  material  (the 
parameters  are  given  on  the  figures ) . 

The  sensitivity  of  resistivity  to  fluence  should  be  noted  near  complete 
compensation.  Furthermore,  in  the  expanded  scale  in  Figure  2,  note  that  the 
maximum  resistivity  occurs  at  a fluence  below  that  necessary  for  100?  compensation. 
This  is  a result  of  the  fact  that  the  carrier  mobilities  are  not  equal. 

By  differentiation  of  Eq.  (13),  it  can  be  shown  that 


2eni^Vp 


where 


is  the  highest  resistivity  obtainable. 

It  can  similarly  be  shown  that  the  fluence,  <t>  , to  obtain  maximum  resistivity, 
is  related  to  the  critical  fluence  to  obtain  100?  compensation  by  the  expression 


"5T 


which  depends  only  on  material  parameters  n^,  the  mobility  ratio,  and  the 
phosphorus  production  rate  K.  This  difference,  therefore,  represents  a 
universal  constant  for  a particular  type  of  material.  It  is  interesting  to  note 
that  since  n^  and  K can  be  measured  independently,  an  experiment  to  generate 
Figure  2 gives  a new  technique  for  measuring  the  drift  mobility  ratios  independent 
of  any  direct  mobility  measurement. 

This  model  can  be  extended  to  include  mixed  conduction  effects  due  to 
compensation  on  the  Hall  factor,  R^,  the  Hall  mobility,  y^,  and  the  apparent 
number  of  free  carriers,  +_  l/R^e. 

From  equations  (ll)  and  (12),  we  found  that 

[l(‘ -s,>Y(‘ l)!  •&)’] 

» - 


Resistivity 


(1  - ♦/♦<■) 

0.277  0.1  0 -0.1 


o * 1.498  X 103  n-cm 
o 

un  = 1500  cro2/v-sec 
up  ■ 600  cm2/v-sec 
♦c  ■ 4.149  X 1016  n/cm2 
«n,  * 4.141  X 1016  n/an2 


n * 1.39  X 10lu  cm'3 


Fluence  (n/cm2) 


6.0  X 10>6 


Figure  1.  Resistivity  and  carrier  concentration 
vs.  fluence. 
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where  (+)  is  for  p>n  and  (-)  is  for  n>p.  This  equation  represents  the  actual 
number  of  carriers  as  a function  of  fluence,  <f>,  and  the  critical  fluence,  <f>c» 
required  for  exact  compensation. 

Because  of  mixed  conduction  effects  near  exact  compensation,  the  Hall 

1* 


2 

np  = n± 

gives 

P(P2  - i 

■A2’ 

e (p2  + 

nA)2 

n(ni2  - 

nV) 

^e(ni2  + 

A)2 

where  b = yn/yp  is  the  ratio  of  drift  mobilities.  Using  the  usual  relation 
between  Hall  mobility  and  R^, 


and  equation  (13) 


^ -_h  ( + ) R>o 

P - 

(-)  v° 


2 2 
eyp(bni  + p ) 


2 2 
eyp(bn  + n1  ) 


y e(bn2  + n.2) 
n l 


we  obtain  0 0 0 

for  Vo 

(p2  + nA) 

yH  = , 2 2 2 « (l8^ 

yn(n  b - ni  ) 

^ b(n2b  + n^2)  and  for 

Note  that  near  complete  compensation,  the  Hall  mobility  is  considerably 
suppressed  as  p + n^b  or  n -*■  n^/b. 

The  apparent  free  carrier  concentration  obtained  from  the  Hall  coefficient 
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(19) 


is,  near  compensation,  larger  than  the  actual  carrier  concentration  since 


PH  = 


(p2  + ni2b)2 


V 


p(p2  - ni2h2) 


°H 


= - 1 
RHe 


(n2!)  + n.2)2 


n(n2b2  - n 2) 


These  equations  are  now  in  a form  from  which  Rjj.  and  or  p^  can  he 
obtained  as  a function  of  fluence  by  direct  substitution  of  Eq.  (l6)  into 
Eqs.  (IT),  (18),  and  (19).  Programs  have  been  written  to  generate  these  variables 
as  a function  of  fluence.  Data  calculated  from  the  code  is  shown  for  a typical 
case  in  Figures  (3),  (U),  (5).  and  (6). 


Note  that  as  p + n^b; 


This  is  not  the  same  condition 


as  obtained  for  the  resistivity  which  reaches  a maximum  at  p 


= bV 


lU 


w 

I 

I 

I 


I 

I 

I 


III.  EXPERIMENTAL  TECHNIQUES  AND  APPARATUS 

In  this  section,  a description  of  the  apparatus  and  techniques  required 
to  irradiate  the  Si  and  to  evaluate  the  material  before  and  after  irradiation 
is  described. 

A.  Reactor  Description 

The  reactor  at  the  University  of  Missouri  is  a high  flux,  light 
water  moderated  research  reactor  designed  to  operate  on  highly  enriched 
fuel.  It  has  a beryllium  reflector,  a graphite  moderator  and  forced 
cooling  system  for  pool  and  core.  The  reactor  is  operated  at  a power  level 
of  10  MW.  The  reactor  runs  routinely  at  10  MW  for  90-100  hours  per  week, 

52  weeks  per  year.  A comparison  of  the  peak  flux  of  the  MURR  reactor  with 
other  research  reactors  in  the  United  States  is  shown  in  Figure  J. 

A cross-sectional  view  of  the  reactor  core  and  various  experimental 
stations  is  shown  in  Figure  8.  One  of  the  major  advantages  of  the  MURR 
reactor  for  Si  irradiation  is  the  easy  access  to  the  sample  positions  in 
the  pool  and  reflector. 

At  present,  two  2"  dia.  x 30"  long  silicon  sample  irradiation  positions 
with  sample  rotators  are  located  in  the  graphite  reflector  in  positions  I- 

and  H-2  in  Figure  8.  The  average  flux  at  these  positions  is  about  3 x 1C* 

? 

n/cm  /sec.  The  pneumatic  tubes  used  for  activation  analysis  rapid 

sampling  handling  are  also  shown  in  Figure  8. 

Beyond  the  graphite  positions  is  the  S-basket  with  seven  3"  dia.  Si 

holes  which  can  be  rotated  manually  or  mechanically.  The  flux  in  the 

13  2 12  2 

S-basket  varies  from  10  n/cm  /sec  to  10  n/cm  /sec. 

For  the  very  low  fluences  required  for  the  exact  compensation  of 

detector  Si,  a variable  flux  pool  facility  was  developed  and  constructed. 

The  samples  in  this  position  are  also  mechanically  rotated.  The  flux 

in  this  facility  is  a maximum  in  the  present  close-in  position  and  is  about 
11  2 

7 x 10  n/cm  /sec.  This  flux  can  be  reduced  considerably  by  moving  the 
facility  further  out  in  the  pool.  The  approximate  position  of  the  bulk 
pool  Si  facility  is  shown  in  Figure  8. 

A typical  flux  profile  for  the  1-2  position  is  shown  in  Figure  9. 
Silicon  ingots  up  to  10"  length  are  positioned  in  the  linear  portion 
of  this  flux  and  flipped  end  for  end  at  precisely  50%  of  the  total  fluence. 
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The  flux  profiles  in  the  other  positions  is  approximately  the  same  shape 
but  tends  to  flatten  at  the  peak  at  positions  further  from  the  core.  In 
the  variable  flux  facility,  the  peak  is  sufficiently  flat  to  make  the  peak 
of  the  flux  profile  useful  over  a 5"  length. 
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B.  Bulk  Pool  Variable  Flux  Irradiation  Facility 


A calculation  of  the  fluence  required  to  compensate  a typical 

12  3 

residual  boron  concentration  of  about  3 x 10  atoms/cm  yields  a fluence 

2 

of  1 to  2 x 10  n/cm  . A typical  loading  or  unloading  time  of  five  minutes 

and  a fluence  accuracy  of  +_  0.5?  restricts  the  irradiation  flux  to  a 

11  12  2 

value  between  1 x 10  to  1 x 10  n/cm  /sec.  Since  this  flux  is  an 
order  of  magnitude  lower  then  that  obtained  in  the  S-Basket,  a special 
irradiation  facility  was  constructed  for  this  project.  The  general 
location  of  this  "Bulk  Pool  Variable  Flux  Facility"  is  shown  in  Figure  8. 

Figure  10  shows  the  construction  details  of  the  moveable  box  which 
is  held  onto  the  S-Basket  by  two  rigid  rods  which  pass  through  holes  A. 

(The  sample  holding  tube  is  horizontal  in  this  diagram  but  vertical  in 
the  pool).  Figure  10  shows  a detail  of  the  support  rods  while  Figure  11 
shows  the  detector  holder  which  inserts  into  the  main  cylinder  in  Figure  10. 
The  sample  holder  is  a cylinder  which  rotates  inside  the  detector  hole 
in  Figure  12.  It  should  be  noted  that  the  neutron  detector  spirals 
concentrically  around  the  sample  position  so  that  a very  accurate 
representation  of  the  flux  on  the  sample  is  continuously  monitored.  The 
thermal  to  fast  flux  in  this  facility  is  estimated  to  be  better  then  30:1. 
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C.  Neutron  Flux  Detection  and  Integration  System 

The  detector  discussed  in  the  last  section  is  a self  powered  detector 
of  the  thin  wire  design. 

103 

Thermal  neutrons  are  captured  in  the  thin  Rh  wire  (100#  abundance). 

A nuclear  reaction  is  induced  which  is  proportional  to  the  flux  of  the 

type  10Vu,  r)  10V . 10kM  . r. 

The  half  life  for  the  0 emission  is  4.2  sec.,  therefore,  this  detector  is 
able  to  follow  the  rather  slow  flux  variations  caused  by  movement  of  the 
reactor  control  rods.  The  electrons  from  this  reaction  are  collected  on 
a concentric  shield  to  produce  a current  proportional  to  the  flux.  This 
current  is  then  integrated  and  an  alarm  set  off  at  a predetermined  preset. 

Considerable  effort  has  been  spent  calibrating  this  system.  Flux 
wires  of  Co  doped  A1  have  been  used  to  measure  the  flux  by  the  radioactivity 
generated  by  a given  fluence  and  compared  to  the  observed  counts  from  the 
integrator  system.  Two  types  of  flux  wire  have  been  used,  N.3.S.  and 
Cohn.  Although  considerable  variation  due  to  flux  wire  counting  statistics 
and  Co  concentration  are  observed,  the  running  mean  of  the  two  sets  of 
flux  wires  agree  to  an  accuracy  approaching  that  of  the  integration  system 
accuracy  ( + 0.2#).  This  is  illustrated  in  Table  4. 

TABLE  4.  Integrator  System  Calibration  Data 
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Wire 

Fluence /Count 
(n/cm2 /count) 

#Devlation  from 
Last ■ Mean 

Running 

Mean 

1 

: 

* 

Cohn 

33-34 

12 

1.557  x 10 

- 0.128# 

1.557  x 10l2 

j 8 

33-34 

1.610 

+ 3.27 

1.594 

9 

33-34 

1.568 

+ 0.577 

1.578 

33- 34B 

1.549 

- 0.641 

1.576 

! 12 

33-34 

1.511 

- 3.079 

1.559 

| 5 

NBS 

12 

1.553  x 10 

- 0.830# 

12 

1.553  x 10X£ 

8 

NBS 

1.482 

- 5.364 

1.518 

11 

NBS 

1.664 

+ 6.258 

1.566 

100  x NBS  - Cohn 

(NBS  + Cohn) /2 

= - 0.449# 
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We  believe  we  have  demonstrated  that  the  fluence  as  determined  by 
the  integration  system  is  accurately  known  and  reproducible  to  about  +0.5%. 


D.  Neutron  Activation  Analysis  System 

Existing  neutron  activation  analysis  (NAA)  facilities  were  employed 

for  measuring  elemental  impurities  in  the  silicon  ingot  starting  materials. 

Two  NAA  irradiation  facilities  are  available:  a pneumatic  tube  transfer 

system  into  a neutron  flux  of  about  8 x 10  n/cm  -sec  which  permits  the 

measurement  of  relatively  short-lived  activation  products;  and  an  outer 

lU  2 

reflector  position  at  approximately  1 x 10  n/cm  -sec  for  long  term 
irradiations.  The  gamma- ray  spectroscopy  instrumentation  consists  of  a 
lead  shielded  50  c.c.  high  resolution  Ge(Li)  detector  coupled  to  a 4096 
channel  MCA  with  computer  compatible  magnetic  tape  output.  An  I.B.M.  370 
computer  is  used  for  data  reduction.  These  facilities  provide  "state  of 
the  art"  capability  and  reliability  for  neutron  activation  analysis  of 
trace  element  impurities  in  silicon. 
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E.  Four-point  and  Thermal -probes 

As  discussed  in  Section  II,  considerable  information  can  be  obtained 
about  the  nature  of  the  compensation  process  by  measuring  resistivity  and 
majority  carrier  type  after  typical  majority  carrier  mobility  is  established 
for  the  ingot  or  wafer. 

Thermal  probe  measurements  of  majority  carrier  type  are  made  in  a 
conventional  manner  using  a variable  temperature  soldering  pencil  tip 
for  the  hot  probe  and  a steel  tip  cold  probe.  The  temperature  gradient 

generated  signal  is  fed  into  a Keithley  Model  602  electrometer'*  operated 

—8 

in  the  center  zero  current  measuring  mode  usually  on  a range  of  10~  to 

_7 

10  A.  The  hot  probe  is  connected  to  the  electrometer  low  while  the 
cold  probe  is  connected  to  the  high  so  that  a positive  deflection  denotes 
p-type  and  a negative  deflection  denotes  n-type.  The  type  of  the  unknown 
sample  is  also  compared  with  two  standard  samples  obtained  from  a multipass 
zone  refined  ingot  (p-type)  and  a neutron  transmuted  and  annealed  sample 
(n-type).  Agreement  of  the  thermal-probe  with  Hall  effect  measurements 
has  been  excellent  up  to  resistivities  approaching  100,000  fi-cm  once 
experience  is  gained  with  the  standard  samples.  Beyond  this  resistivity, 
mixed  typing  is  often  encountered  suggesting  that  nonuniform  doping  is 
observed  before  the  disappearance  of  the  thermoprobe  signal  is  observed. 

Several  four  point  probe  heads  of  different  probe  specing  and  probe 
pressure  have  been  tried  with  excellent  results  in  resistivity  ranges 
usually  inaccessible  to  four-point  probe  measurement.  The  system  consists 
of  a Keithley  Model  225  constant  current  source'* , a Hewlett  Packard  Model 
3I+65A  digital  multimeter  , and  an  Alessi  Industries  ATP-1  Test  Probe  Fixture  1 
which  has  been  modified  to  have  a 12"  vertical  travel  and  a sample  translation 
stage.  Two  types  of  four-point  probe  heads  have  been  tried,  an  Alessi 
Model  Al+P-25  with  a probe  spacing  of  0.025",  a probe  load  of  70-180  gm;  and 

g 

a Signatone  probe  with  a spacing  of  O.O65"  a tip  radius  of  0.1  mil  and  a 
probe  load  of  180  gms.  Both  probe  heads  have  tungsten  carbide  tips.  Although 
both  probe  heads  work  well,  the  larger  probe  spacing  is  more  forgiving  in 
current  selection  and  care  of  measurement  as  has  been  used  except  for  very 
small  samples  where  edge  effects  become  a problem.  The  usual  probe  correction 
factors  have  been  applied  to  obtain  the  resistivity  from  the  measured 
resistance. ^ ^ 
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To  account  for  nonunifonr.  probe  spacing,  the  probes  were  lowered  onto 

an  aluminized  microscope  slide  and  the  spacing  measured  with  a traveling 

microscope.  This  procedure  was  repeated  several  times  and  a correction 

9 

term  generated  from  theory.  This  correction  factor  was  also  checked 

by  measuring  the  resistivity  of  a transmutation  doped  and  annealed  sample 

12 

which  had  been  calibrated  by  the  Van  der  Paw  technique.  Agreement  was 
found  to  be  better  than  1%  for  a 5 ft-cm  sample.  All  samples  are  measured 
in  a dark  box. 

> 

Since  the  A.S.T.M.  Standard  for  resistivities  over  2000  ft-cm  suggested 
that  measurements  could  not  be  made  to  an  accuracy  of  better  than  +15$  , 
extreme  care  was  taken  to  check  the  probe  calibration  factor  periodically 
for  high  resistivity  measurements. 

An  ingot  which  had  been  cut  into  five  bulk  samples  (p  'u  2500  fd-cm) 
and  characterized  by  four  point  probes  at  Rockwell  International  was 
measured  at  MURR.  Care  was  taken  to  use  a current  sufficiently  low  that 
sample  heating  did  not  affect  the  measurements.  The  results  of  I-V  plots 
for  two  sides  of  sample  VZ-072-#5  are  shown  in  Figures  13  and  1*4.  The 
nonlinearity  in  these  measurements  at  high  currents  is  due  to  sample 
heating. 

Note  that  the  linear  portions  of  the  curves  do  not  go  through  the 
origin.  This  offset  voltage  can  be  measured  on  the  voltmeter  in  the  absence 
of  current  and  is  observed  to  vary  with  probe  pressure,  surface  condition 
of  the  sample  and  position  of  the  probes  on  the  sample.  Although  this 
effect  is  negligible  for  5 ft-cm  material,  the  offset  voltage  can  become  an 
appreciable  part  of  the  signal  voltage  for  high  resistivity  samples.  It 
was  also  noticed  that  the  slope  near  the  origin  is  independent  of  the  offset. 

This  fact  allows  us  to  make  correct  measurements  up  to  very  high  resistivities 
if  proper  data  handling  is  used.  The  problem  is  overcome  by  measuring  both 
the  positive  and  negative  voltages  resulting  from  positive  and  negative 
supply  currents  to  determine  the  effective  resistance.  This  resistance  is 
then  equal  to  [ V(+)  - V(-)  ] / [l(+)  - I(-)  ].  This  resistance  determines 
the  slope  of  the  lines  in  Figures  (13 ) and  (lM  and  the  product  of  this 
slope  times  the  probe  factor  determines  the  resistivity. 

Results  of  a series  of  measurements  on  the  five  bulk  samples,  VZ-072-01 
through  # 5,  are  compared  to  results  obtained  by  Dr.  Arst  at  Rockwell  and 
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ar**  shewn  in  Table  5. 


TABLE  5*  A comparsion  of  four  point  probe  resistivity  measurements  made 
at  Rockwell  and  MURR. 


Sample 

Rockwell  Data 
Resistivity  (fi-cm) 

MURR  Data 
Resistivity  (fi-cm) 

100  Pr-Pm 
Pra 

0l(T) 

2.29  x 103 

2.27  x 102 

+0.88 

01(B) 

2.12 

2.1*3 

-0.1*1 

02(T) 

2.1*3 

2.1+5 

-0.82 

02(B) 

2.53 

2.53 

0.00 

#3(T) 

2.50 

2.52 

-0.99 

03(B) 

2.57 

2.58 

-0.39 

0U(T) 

2.57 

2.58 

-0.39 

0MB) 

2.55 

2.56 

-0.39 

05  (T) 

2.55 

2.55 

0.00 

1 

05(B) 

2.53 

2.1*9 

+1.61 

It  is  seen  that  the  agreement  between  the  two  laboratories  is  quite  good 
for  this  resistivity  range  the  average  deviation  being  about  +_  0.6%.  The 
agreement  between  the  four  point  probe  measurements  has  been  observed  to 
decrease  with  further  increases  in  resistivity.  On  comparing  both  laboratories 
four  point  probe  measurements  with  Van  der  Pauw  measurements  on  samples  with 
resistivities  in  the  7000-8000  fi-cm  range  we  find  the  deviation  at  MURR 


from  Van  der  Pauw  to  be  +_  1.5%  while  the  deviation  of  the  Rockwell  measurements 
j | is  + 5%. 

This  increased  error  with  resistivity  seems  to  be  following  a power 


law  dependence  given  by  the  expression 


100  — 
p 

■ = 0.1  p0-71 

(A.S.T.M.) 

100  - 
p 

■ = 0.001  p0'81 

(MURR) . 

The  physical  meaning  of  these  equations  is  obscure,  however,  a 
similar  power  law  is  observed  for  the  calibration  curves  for  spreading 
resistance  measurements.  The  difference  between  the  A.S.T.M.  and  MURR 
values  for  K are  not  known,  however,  the  A.S.T.M.  data  is  taken  for  both 
n-  and  p-type  Si  in  various  orientations  while  the  MURR  data  is  taken 
on  p-type  only  on  a (ill)  face.  The  upper  expression  is  useful  for 
estimating  the  maximum  expected  error  for  intrinsic  material  at  230,000 
fi-cm.  We  find  an  expected  value  of  about  +_  22$  at  this  resistivity  value. 

This  estimated  error  is  in  reasonable  agreement  with  deviations  in  resistivity 
for  samples  neutron  irradiated  to  the  same  fluence  but  before  annealing 
(nearly  intrinsic  but  slightly  p-type). 
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Minority  Carrier  Lifetime  Apparatus 
The  minority  carrier  lifetime  of  all  ingots  before  irradiation  and 
some  after  irradiation  and  annealing  have  been  measured  by  the  photoconductive 
decay  technique.  Ingots  are  coated  with  Aquadag  solution  on  each  end  and  held 
in  a small  lab  vise  to  which  braided  copper  shielding  wire  has  been  flattened 
and  attached  to  the  vise  Jaws  for  contacting.  The  contacting  area  is 
shielded  from  the  light  to  prevent  photovoltaic  effects  from  the  contacts. 

This  effect  can  be  detected  if  observed  as  a linear  decay  of  the  Junction 
capacitance  of  a blocking  contact. 

A small  power  supply  supplies  the  bias  voltage  across  the  sample  and 
a variable  series  load  resistor.  This  resistor  is  adjusted  for  maximum 


signal  before  data  is  taken,  an  essential  procedure  for  interpretation  of 
It 

the  decay.  Optical  excitation  is  provided  by  a General  Radio  Strobotac 
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model  1531-AB  whose  pulse  time  is  short  compared  to  the  photoconductive 


decay  time  constant.  The  signal  is  observed  with  a Tektronix  Type  51*? 
oscilloscope  with  a dual  trace  plug-in  unit.^ 


Typical  decay  voltage  as  a function  of  time  are  shown  in  Figure  15. 


It  can  be  shown  that  the  voltage  signal  as  a function  of  time  for  the 

IT 


condition  that  R load  = R sample  is  given  by 


V(t)  = 


-t/x 


Vn  / 

' © ' o 


(20) 


bias 


1 + [i  it) 


-t/it 


where  (An/n°)0  is  the  maximum  change  in  minority  carrier  concentration 


and  x is  the  carrier  lifetime.  For  times  long  compared  to  the  optical 
excitation  decay 


V(t) 


V, 


(ft 


-t/x 


Therefore,  a semilog  plot  of  V(t)/V  has  a slope  of  x 

D X £LS 


bias  ' o ' o 

where  x is  the 

U-LO.& 

minority  carrier  lifetime. 

Care  has  been  taken  not  to  optically  overmodulate  the  samples  since 
this  is  known  to  produce  exaggerated  lifetimes.  The  procedure  we  use  is 
to  reduce  the  light  intensity  until  the  observed  decay  time  constant  becomes 
a minimum.  It  is  also  important  that  the  carrier  sweep  out  time  given 


UV.  . 
bias 
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where  L is  the  sample  length  and  p the  carrier  mobility,  be  long  compared 
to  the  decay  time. 

Trapping  effects  have  been  detected  in  some  samples  by  observing  a 
decrease  in  the  decay  time  when  the  sample  is  simultaneously  illuminated  from 

a dc  tungsten  light  source  (4m><E  ).  In  these  cases,  the  lifetime  is  assumed 

& 

to  be  that  observed  during  dc  illumination  if  the  reduction  in  decay  time 
saturates. 

Circuitry  has  been  developed  to  eliminate  the  need  to  take  data  from 
the  exponential  photoconductive  decay  displayed  on  the  oscilloscope.  A 
variable  time  constant,  variable  pulse  height  exponential  generator  has 
been  constructed  which  provides  a decay  time  constant  whose  value  can  be 
obtained  from  the  setting  of  a ten-turn  pot.  This  exponential  pulse  is 
triggered  from  the  strobotac,  inverted,  and  added  to  the  photoconductive 
decay  signal.  The  oscilloscope  then  functions  as  a null  detector.  A 
straight  line  signal  indicates  that  the  exponential  generator  has  matched 
the  photoconductive  signal  in  both  amplitude  and  time  constant.  Deviations 
from  non-exponential  decays  and  trapping  effects  are  easily  detected  with 
this  system.  A schematic  of  the  system  is  shown  in  Figure  16  while  a 
circuit  diagram  for  the  exponential  pulse  generator  is  shown  in  Figure  IT. 


r 

i . 


i } 

t j 


39 


CARRIER  LIFETIME  MEASUREMENT  APPARATUS 


EXPONENTIAL  DECAY  GENERATOR 


C-.  High  Impedance  Van  der  Pauw  Resistivity  and  Hall  System 
(Electronics,  Vacuum  System,  and  Helitran) . 

Resistivity  and  mobility  can  be  determined  as  a function  of  temperature 

from  about  8°K  to  room  temperature  by  this  system.  Temperature  control  is 

1 8 

obtained  from  a Helitran  and  temperature  controller.  The  system  is 

capable  of  measuring  resistances  of  high  impedance  samples  up  to  a 

12  19 

resistivity  of  about  10  fi-cm.  A dc  guarded  lead  approach  in  which 

each  sample  lead  guard  is  driven  by  a separate  high  input  impedance  unity 

gain  amplifier  (Keithley'’  Model  602)  to  the  voltage  on  the  inner  lead  reduces 

leakage  currents  and  improves  the  system's  time  constant.  The  circuitry 

for  this  approach  is  shown  in  Reference  19.  The  digital  voltmeter  and 

constant  current  source  from  the  four  point  probe  system  are  also  used  in 

this  system.  A switching  circuit  connecting  the  appropriate  sample  leads 

to  the  proper  electrometer  and  current  source  allows  either  conventional 

12 

resistivity  and  Hall  effect  or  Van  der  Pauw  measurements  to  be  made. 

This  circuit  is  shown  in  Figure  18. 

A vacuum  system  has  been  constructed  to  provide  evacuation  for  the 

Helitran  system  and  for  the  vacuum  annealing.  The  system  uses  an  oil 

diffusion  pump,  mechanical  forepump  and  liquid  nitrogen  cold  trap.  This 

_7 

system  provides  a vacuum  of  5 x 10  torr. 

The  magnetic  field  for  Hall  effect  measurements  is  provided  by  a Varian 

20 

Model  VljOOk  four-inch  magnet  system.  The  magnetic  field  has  been 
calibrated  to  + 1.3%  using  both  a rotating  coil  gauss  meter  which  has  been 
compared  to  an  NMR  probe  and  a calibrated  semiconductor  Hall  probe.  The 
semiconductor  Hall  probe  is  used  to  set  the  magnetic  field  using  circuitry 
which  provides  digital  readout  of  the  field.  This  circuitry  is  shown  in 
Figures  19  and  20.  The  magnetic  field  for  Hall  effect  measurements  in 
this  report  is  6.195  kG. 

Although  it  would  be  desirable  to  determine  the  dominate  energy  levels 
at  various  points  along  the  annealing  curves  by  a measurement  of  n vs  T, 
we  have  not  been  able  to  obtain  this  data  yet  due  to  the  very  poor  customer 
service  provided  by  Air  Products.  Although  the  cryogenic  part  of  the 
Helitran  system  has  functioned  properly  from  the  time  it  was  received, 
the  whole  system  was  delivered  six  weeks  after  the  date  promised.  Furthermore, 
another  six  week  delay  resulted  from  improper  Helitran  wiring  which  required 


12 


DIGITAL  KGAUSS  METER 


1 


Ik 


a return  to  the  manufacturer.  An  additional  three  month  delay  was  caused 
by  late  delivery  of  the  temperature  controller  and  its  failure  on  the 
second  cool-down.  Finally,  a heater  short  to  a thermalcouple  was  discovered 
which  caused  heater  failure  and  another  two  week  delay;  a total  delay  of 
nearly  five  months  beyond  the  delivery  date  promised  by  Air  Products  in 
their  bid  package. 

We  can  only  apologize  for  this  situation  over  which  we  have  had  little 
control  and  report  that  we  have  made  every  effort  to  get  the  attention  of 
the  management  of  Air  Products  (over  15  phone  conversations  and  a series 
of  letters  to  every  level  of  the  corporation)  with  little  success  except 
for  a phone  call,  by  them,  to  ask  when  we  intend  to  pay  for  the  Helitran! 

We  believe  the  system  to  be  functioning  properly  at  this  time. 
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H.  Photoconductivity  Syster. 

Because  the  transmutation  doping  studies  have  been  performed  thus 
far  on  undoped  material,  the  necessity  of  performing  photoconductivity 
studies  to  obtain  the  degree  of  compensation  or  the  photoconductive  spectral 
response  have  not  been  needed.  This  system  has  been  developed  in  anticipation 
of  work  hopefully  to  be  completed  near  the  end  of  the  contract.  The  system 
consists  of  a glow  bar  light  source,  a Perkin-Elmer  double  pass  prism 
monochrometer  and  electronics,  two  thermocouple  detectors  and  exit  optics. 

The  prism  and  windows  are  KRS-5  material  ( 1—50  Urn).  The  exit  beam  from 
the  monochrometer  which  is  diverging  is  intercepted  by  a spherical  mirror 
to  produce  a parallel  beam  of  light.  Part  of  this  parallel  beam  is  interrupted 
by  a diagonal  plane  mirror  and  focused  on  a reference  thermocouple  detector 
with  a KRS-5  window.  The  other  part  of  the  beam  is  focused  through  a 
second  window  onto  the  photodetector  sample  mounted  in  the  Helitran  by 
another  spherical  mirror.  By  taking  ratios  of  these  two  signals  the 
spectral  response  of  the  photodetector  can  be  obtained  normalized  to  the 
light  source  and  system  response.  The  final  construction  of  the  output 
optics  table  is  now  underway  and  should  be  completed  soon. 


I . Miscellaneous:  Sample  Preparation  Equipment,  Annealing  System, 

Irradiation  of  Silicon,  Post-Irradiation  Handling 

Samples  are  cut  into  wafers  either  by  a diamond  wire  saw  or  diamond 
disk,  cut-off  saw.  The  latter  has  been  found  more  suitable  for  silicon. 

Samples  are  heavily  etched  after  cutting  to  remove  surface  damage  in  a 
conventional  HF/HNO^/acetic  etch.  Samples  are  then  ground  using  various 
grades  of  SiC  and  Al^O^  grits  and  finally  washed  in  deionized  water, 
acetone,  methanol,  and  trichloroethylene. 

Samples  are  annealed  in  a tube  furnace  inside  a Spectrasil  quartz 

_r7 

tube.  The  sample  may  be  annealed  in  either  a vacuum  of  5 x 10  ‘torr.  or 
in  a flowing  inert  gas  atmosphere.  Sample  surfaces  are  protected  on  each 
side  with  a wafer  of  float  zone  silicon  to  prevent  contamination  by  the 
quartz.  The  furnace  is  capable  of  temperatures  of  up  to  1000°C.  The 
temperature  is  regulated  to  +_  0.5°C. 

Samples  to  be  irradiated  are  placed  between  silicon  spacers;  the  whole 
stack  of  samples  and  spacers  is  then  wrapped  in  aluminum  foil.  The  flux 
wires  are  now  wrapped  around  the  foil  and  the  whole  package  is  placed  in 
an  aluminum  can  which  is  welded  shut  to  protect  the  silicon  from  the  pool 
water.  Radioactivity  is  monitored  after  a cool-off  period  in  the  reactor 
pool.  For  the  very  light  irradiations  in  this  project,  radioactivity  decays 
to  background  very  quickly  and  has  usually  been  attributed  to  surface 
contamination  where  encountered. 

A considerable  effort  has  been  spent  studying  surface  contamination 
after  irradiation  but  before  annealing.  By  careful  handling,  no  surface 
contamination  has  been  detected  before  annealing  except  for  a Au  contamination 
on  one  sample  from  sawing  in  another  laboratory  before  our  saws  were 
operational.  A few  samples  have  also  -shown  traces  of  Na  from  handling. 

This  has  been  eliminated  by  a light  HF  cleaning  after  irradiation  but 
before  annealing.  We  are  confident  that  our  samples  are  extremenly  clean 
prior  to  annealing  based  on  activation  analysis. 
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IV.  EXPERIMENTAL  RESULTS  AMD  DISCUSSION 

In  this  section,  the  prior-to-irradiation  material  evaluation,  irradiation 
and  annealing  experiments  are  discussed.  We  have  noted  several  differences  in 
our  resuits  from  previously  published  work  and  will  discuss  these  differences 
in  the  following  subsections.  A preliminary  activation  energy  for  the  electrical 
activity  of  the  irradiation  produced  donors  has  been  measured.  We  have  been 
able  to  determine  that  the  best  single  annealing  temperature  for  NTD-Si  at 
above  800°C  for  one  hour  followed  by  a very  slow  cool  down. 

Although  material  from  a variety  of  Si  producers  has  been  evaluated,  we 
have  found  that  control  by  compensation  on  the  p-type  side  of  type  conversion 
from  5000  ft-cm  to  intrinsic  is  impossible  in  the  material  we  have  investigated. 

We  will  present  data  which  suggests  that  residual  oxygen  is  the  source  of 
difficulty. 

We  have  also  demonstrated  very  good  control  of  resistivity  on  the  n-type 
side  of  type  conversion.  We  have  shown  that  the  final  resistivity  can  be 
controlled  to  at  least  +2%  using  only  a preliminary  check  of  ingot  mobility  and 
four  point  probe  measurements.  We  have  sound  experimental  evidence  at  this 
point  which  suggests  that  it  may  be  possible  to  achieve  this  degree  of  control 

up  to  n-type  resistivities  as  high  as  50,000  to  70,000  ft-cm.  This  corresponds 

—2  12  3 
to  an  excess  donor  concentration  of  ^2  x 10  ppb  or  about  1 x 10  ^ donors/cm“. 

It  may  thus  be  possible  to  compensate  residual  boron  concentrations  of  about  this 

order  of  magnitude  or  larger. 

One  of  the  more  exciting  possibilities  to  develop  out  of  the  present  work 
is  residual  oxygen  decoration  by  damage  defects.  Circumstantial  evideree  suggests 
that  oxygen  detection  in  Si  can  be  pushed  at  least  three  orders  of  magni  ude 
below  present  detection  limits.  If  the  preliminary  results  can  be  substantiated 
by  EPF.  and  optical  measurements,  considerable  progress  in  Si  purification  can  be 
made  as  a result  of  this  technique. 
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A. 


Materials  Evaluation  Before  Irradiation 


All  the  material  investigated  has  been  multizone  refined  float  zone 
silicon  with  resistivities  from  2500  Q-cm  up  to  25000  fi-cn.  The  material 
is  p-type  and  not  intentionally  doped.  The  dominant  acceptor  is  un- 
do ubt  ably  boron.  All  of  the  material  is  nominally  dislocation  free  and 
pulled  in  a (ill)  direction  except  the  Rockwell  material  which  is  dislocated 
and  pulled  in  a (100)  direction.  Neutron  reflection  rocking  curves  show 
reflection  line  widths  of  about  1 sec  of  arc  for  the  dislocation  free 
material  and  about  1 min  of  arc  for  the  dislocated  material.  Material 
investigated  has  come  from  Dr.  Arst,  Rockwell;  Topsil,  Wackers,  Monsanto, 
and  a Texas  Instruments  Lopex  sample  supplied  by  Dr.  Spry  of  AFML.  We 
have  added  an  in-house  number  to  each  ingot  or  wafer  (MURR  #)  and  have 
cross-referenced  these  numbers  with  whatever  sample  designations  have 
been  assigned  to  the  samples  by  the  suppliers. 

Much  effort  has  been  expended  in  trying  to  purchase  float  zone  silicon 

which  is  doped  with  either  In  or  Ga.  Approximately  thirty  manufacturers 

of  silicon  were  contacted  by  phone  and  requests  for  bids  sent  out  to  ten, 

however,  no  bids  were  received  for  this  type  of  material.  We  have,  however, 

received  free,  doped  samples  from  a number  of  sources  and  greatfully 

acknowledge  their  help.  We  have  received  Czochralski  grown  Si:Ga  doped 

l6  1 8 3 

to  a variety  of  levels  in  the  range  of  10  to  10  cm-J  from  Drs.  M. 
Bressler  and  C.  M.  Parry,  Aerojet  ElectroSystems  Co.  We  have  also  re- 
ceived a supply  of  Wacker's  Ga  doped  vacuum  float  zone  material  from  Dr. 

M.  Schlacter,  Aeronutronic  Ford.  Dr.  Arst  of  Rockwell  has  also  supplied 
us  with  Ga  doped  float  zone. 

Because  of  difficulties  in  obtaining  doped  material  and  because  of 
the  poor  service  supplied  to  us  by  Air  Products,  we  have  not  as  yet 
characterized  the  Ga  doped  material  sufficiently  to  irradiate  it,  however, 
this  phase  of  the  program  has  now  begun  and  hopefully  will  be  completed 
by  the  end  of  the  contract. 

Table  6 summarizes  sample  characterization  before  irradiation.  As 
discussed  in  a previous  section,  four  point  resistivity  measurements  have 
been  compared  with  Van  der  Pauw  resistivity  measurements  for  calibration 
purposes.  Mobility  measurements  have  been  taken  from  one  slice  of  each 
ingot  and  the  values  are  very  near  average  values  giver,  in  Table  3 • 
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TABLE  6.  Summary  of  Undoped  Sample  Characterization  Before  Irradiation 


MURR  n 

Source 

P(l+  pt) 
(ft-cm) 

P (V.d.P. 
(ft- cm) 

) 

(cm2/V-sec) 

Type 

T.P. 

Hall 

T 

( psec) 

0 

Monsanto  FZ 

7^9.0 

P 

- 

1 

Rockwell 

VZ-072-3 

2560. 

P 

110 

1 2 

VZ-072-3 

2953. 

P 

110 

3 

VZ-072-3 

3156. 

P 

110 

* 

VZ-072-3 

2721. 

1+61+ . 75 

P 

P 

110 

I 5 

Monsanto  FZ 

1690. 

1383.1+ 

N 

N 

- 

6 

Monsanto  FZ 

1037.5 

ll+OU.5 

N 

N 

- 

i 7 

Rockwell 

VZ-072-2 

21+50- 

2530  (ingot) 

P 

112 

i 

Rockwell 

VZ-083-1+ 

7630. 

7729. 

521+.8 

P 

P 

_ 

9 

| 

Rockwell 

VZ-083-5 

8189. 

8067. 

527.8 

P 

_ 

* 10 

Monsanto  CZ 

(after  irrad 
& anneal) 

3.222 

1267.0 

S 

N 

- 

1 

Monsanto  FZ 

1376.2 

11+73. 5 

N 

N 

- 

Ih 

1 

Topsil 

36-520-0-111 

P 

1200 

15 

Topsil 

36-520-0-111 

1+710- 

6167 

(l"  ingot) 

P 

1200 

| 1 6 

1 

Topsil 

36-520-0-111 

1+71+2. 

P 

1200 

17 

{ 

* 

Wackers 

30661-6 

66512.  ? 

P 

530 

18 

1 

T.I.  Lopex 
LI  09  2 

1+607 

P 

250 

j 19 

Topsil 

36-52O-O-III 

I+9I+2 . 

1+60. 

P 

P 

1200 

20 

T.I.  Lopex 
L1092 

5195.9 

1+1+3.86 

P 

P 

250 

1 

Wacker 

30661-6 

19370. 

P 

530 

(trapping) 

j 22 

Wacker 

30661-6 

201+20. 

1+27.0 

P 

P 

530 

J"  11111 11  111  111  W I'1"11"- 

TABLE  6.  Summary  of  Undoped  Sample  Characterization  Before  Irradiation  (Cont.) 


MURR  P 

Source 

P (1+  pt) 
(ft-cm) 

Kv.d.p. : 

(ft-cm) 

1 WD 

( cm2 /V-sec ) 

Type 

T.P.  Hall 

T 

( ysec) 

23 

Similar  to  MURR 

10 

N 

- 

2k 

Wacker 

30661-6 

For  NM 

P 

530 

25 

T.I.  Lopex 
L1092 

For  NAA 

P 

250 

26 

Topsil 

36-52O-O-III 

For  NAA 

P 

1200 

27 

Topsil 

36-52O-O-III 

52U8 

P 

1200 

23 

Monsanto  FZ 

15^0 

N 

- 

29 

Ro ckwell 
VZ-083-17 

12509. 

P 

_ 

30 

Rockwell 

VZ-083-18 

12796 

F 

_ 

31 

Rockwell 

VZ-083-19 

13326 

P 

_ 

32 

Rockwell 

VZ-083-20 

13687 

P 

_ 

33 

Rockwell 

VZ-083-21 

ll+l60 

P 

. 

3U 

Rockwell 

VZ-083-22 

11+382 

P 

__ 

35 

T.I.  Lopex 
L1092 

5895 

P 

36 

T.I.  Lopex 
L1092 

5587.6 

P 

- 

To 

r VZ-072-1 

2270-21+30 

Check 

Uniformity  * 

VZ-072-2 
' . VZ-072-3 
VZ-072-U 

21+50-2530 

2520-2580 

2580-2560 

VZ-072-5 

2550-21+90 

v 
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Note,  however,  the  Wacker's  sample  is  some  13$  too  low.  This  is 
probably  due  to  the  effect  of  mixed  conduction  in  this  very  high  resistivity 
sample. 

Figure  21  shows  the  axial  resistivity  and  minority  carrier  lifetime 
of  ingot  VZ-072  after  cutting  it  into  slices  of  approximately  1"  thickness. 
Figures  22  through  2?  show  the  radial  uniformity  of  these  samples  as  obtained 
from  four  point  measurements. 

9 

The  rise  in  resistivity  is  due  to  edge  effects  and  is  not  real. 

Table  7 summarizes  the  minority  carrier  lifetime  for  the  various 
ingots  which  have  been  measured. 

TABLE  7.  Summary  of  Minority  Carrier  Lifetimes 


Source 

Ingot  No. 

T Measured 

T Specified  by  Manuf, 

Rockwell 

VZ-072 

100  ysec 

100  ysec 

Topsil 

36-52O-O-III 

1200  (trapping) 

1700 

T.I.  Lopex 

L 1092 

250 

500 

Wacker 

30661-6 

530 

3200 

We  believe  that  the  lifetimes  in  this  resistivity  range  are  probably  over 
exaggerated  by  the  suppliers  because  of  overmodulation  in  experimental 
set  ups  designed  in  general  for  lower  resistivity  material.  It  should  be 
noted  that  our  observed  lifetimes  for  the  Rockwell  material  are  ir  general 
agreement  with  their  measurements.  The  lifetime  of  the  Topsil  material 
indicates  some  evidence  of  trapping  which  does  not  saturate.  It  is, 
therefore,  possible  that  the  real  lifetime  is  shorter.  We  would  estimate, 
however,  that  it  is  probably  not  shorter  than  500ysec. 

Table  8 summarizes  the  drift  mobilities  observed  for  the  various 
ingots  obtained  from  Hall  mobilities  using  the  relationship  for  the  Hall 
factor  shown  in  the  table. 
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TABLE  8.  Summary  of  Drift  Mobilities 


Source 

Ingot  No. 

Type 

Up  (cm2/V-sec) 

Rockwell 

VZ-072 

P 

U65 

Rockwell 

VZ-082 

P 

525-528 

Topsil 

36-520-0-111 

P 

1+60 

T.  I.  Lopex 

L1092 

P 

U1*U 

Wacker 

30661-6 

P 

U27 

Monsanto  FZ 

N 

1383-1*473 

Typical  from  Table 

3 

P 

U9L 

N 

1396 

rP  = YyD  = 

0.77 

rN  = X-3 

We  conclude,  based  on  the  data  in  Tables  7 and  8 that  the  materials  from 
these  sources  is  of  comparable  crystal  quality.  The  Rockwell  material  is  slightly 
inferior,  however,  the  material  was  pulled  in  a (100)  direction  which  is  more 
difficult.  Apparently,  a (ill)  Rockwell  pull  would  be  of  comparable  crystal 
quality.  In  addition,  the  Rockwell  material  has  the  added  advantage  of  extremely 
uniform  boron  distribution  as  evidenced  by  Figures  22-27. 


01  stance  of  Probe  Center  From  Left  Edge  of  Sample  (In) 


Dlsta 


Figure  27.  Radial  resistivity  of  VZ-072-5-bottom 
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B.  Neutron  Activation  Analysis 

An  instrumental  neutron  activation  analysis  (INAA)  procedure  for 

measuring  trace  element  imparities  in  high  purity  silicon  was  established. 

This  INAA  procedure  serves  as  a screening  technique  for  a wide  variety  of 

unexpected  impurity  elements  and  also  provides  quantitative  upper  limit 

data  for  a number  of  important  impurity  elements  if  they  are  not  observed. 

This  latter  group  of  important  elements  include  (a)  elements  known  to 

produce  striking  effects  on  the  physico-chemical  properties  of  silicon, 

(b)  elements  which  are  known  to  be  common  impurities  in  technical  grade  < 

silicon,  and  (c)  elements  which  may  serve  as  indicators  for  contamination 

of  the  silicon  ingots  during  refining  and  handling. 

Short  Irradiation  Procedure.  A thick  wafer  weighing  approximately 

1 gram  was  etched  with  HF-HNO^  and  rinsed  with  deionized  water  to  remove 

any  surface  contamination.  The  cleaned  sample  was  wrapped  in  a precleaned 

13 

polyethylene  bag,  and  irradiated  for  20  min.  at  a neutron  flux  of  C x 10 

n/cm^-sec  using  the  pneumatic  tube  transfer  facility.  Following  irradiation 

the  sample  was  washed  with  hot  8N  HNO^  and  rinsed  with  deionized  water. 

The  gamma-ray  spectrum  of  the  irradiated  sample  was  recorded  for  2000 

seconds  after  a decay  period  of  60  minutes  using  a Ge(Li)  detector.  A 

3 millimeter  Fb  absorber  was  placed  between  the  sample  and  detector  to 
31 

reduce  the  Si  beta  bremstrahlung.  Elemental  standards  prepared  from 
high  purity  elements  or  weighable  compounds  were  irradiated  and  counted 
in  the  same  manner  as  the  samples. 

Long  Irradiation  Procedure.  Following  the  short  irradiation  procedure, 

the  same  samples  were  wrapped  in  aluminum  foil,  placed  in  a reflector 

irradiation  canister,  and  irradiated  for  1+0  hours  at  a neutron  flux  of 

l4  2 

approximately  1 x 10  n/cm  -sec.  After  irradiation  the  samples  were  etched 
with  HF-HNC^  and  rinsed  with  deionized  water  to  remove  surface  contamination 
from  the  aluminum  foil  wrapping.  Each  sample  was  counted  for  2000  sec  after 
decay  periods  of  approximately  ko  hours  and  5 days.  The  gamma-ray  spectra  fj 

from  0 to  h MeV  was  recorded  using  a Ge(Li)  detector.  Elemental  standards 
were  irradiated  and  counted  according  to  the  same  schedule.  The  analytical 
conditions  are  summarized  in  Table  1. 

Data  Reduction.  Existing  computer  codes  were  used  to  reduce  the 
gamma-ray  spectral  data  for  both  samples  and  standards  to  peak  energies 
and  peak  areas  for  all  observed  gamma- rays.  Isotope  assignments  were 


based  on  calculated  energies  and.  where  necessary,  half-life  measurements . 
Elemental  concentrations  for  observed  impurities  were  calculated  by 
comparison  of  sample  and  standard  peak  areas  after  proper  normalization  of 
differences  in  irradiation  exposure  and  decay  time,  if  necessary.  Minimum 
detectable  peak  areas  for  isotopes  not  observed  were  estimated  as  three 
times  the  statistical  deviation  of  the  spectral  base  line  in  the  peak 
region.  Detection  limits  for  the  elements  measured  are  given  in  Table  9 
for  a 1 gram  silicon  sample  using  a value  of  2.3^  gm/c.c.  for  the  density 
of  silicon. 

Results . Silicon  ingots  obtained  from  four  different  suppliers  were 

analyzed.  Elemental  concentration  values  are  given  in  Table  10  as  atoms 

of  impurity  per  c.c.  of  silicon  assuming  that  the  impurity  is  homogeneous 

in  the  matrix.  The  values  given  for  sodium  probably  represent  surface 

27  2U 

contamination  derived  from  the  A1  (n,  a)  Na  reaction  occurring  in  the 
aluminum  foil  encapsulation. 

Figure  28  illustrates  a sawing  contamination  which  was  detected  on 
VZ-072.  This  sample  was  cut  on  a glass  cut-off  tool  outside  our  laboratories. 
It  is  apparent  that  gold  surface  contamination  was  present  and  was 
difficult  to  remove  by  etching.  The  etching  procedure  apparently  tends 
to  selectively  etch  around  the  gold  spots  and  undercut  until  the  gold 
is  removed. 
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TABLE  10.  Elemental  Concentrations  Measured  in  Four  Silicon  Ingots 
Impurity  level  measured  in  atoms  of  impurity/c . c . * 


.enent 

Rockwell  VZ-072 

Topsil  36-520 

Wackers  30661-6 

TI  Lopex  L1092 

In 

n.d. 

n.d. 

n.d. 

n.  d. 

Mn 

n.  d. 

12 

^2  x 10 

n.d. 

n.  d. 

Mo 

1 P 

,'J*  x 10 

n.d. 

n.  d. 

n.d. 

Cu 

12 

2.9  x 10 

Vf  x 1011 

12 

1.1  x 10 

12 

1.3  x 10 

Ga 

n.d. 

n.d. 

n.d. 

n.d. 

Au 

4.7  x 109 

n.d. 

5.3  x 109 

9 

^3  x 10* 

W 

n.  d. 

n.d. 

n.d. 

n.d. 

As 

n.d. 

n.d. 

n.d. 

n.d. 

Sb 

n.d. 

n.d. 

n.d. 

n.d. 

Yb 

n.d. 

n.  d. 

n.d. 

n.d. 

Na 

8.9  x 1013 

13 

2.0  x 10  J 

1.8  x 1013 

13 

1.9  x 10  ^ 

Eu 

n.d. 

n.  d. 

n.  d. 

n.  d. 

Cr 

n.  d. 

n.d. 

n.d. 

n.  d. 

Th 

n.d. 

n.  d. 

n.d. 

n.  a. 

Sc 

n.  d. 

n.  d. 

n.  d. 

n.d. 

Hf 

n.d. 

n.d. 

n.  d. 

n.d. 

Co 

n.  d. 

n.  d. 

n.d. 

n.d. 

Ag 

n.d. 

n.  d. 

n.d. 

n.d. 

Ta 

n.d. 

n.d. 

n.d. 

n.d. 

* n.d. 

indicates  not  detected 

concentrations 

above  the  background.  Detection 

limits  are  given  in  Table  9- 


m ^ 


C.  Irradiation  and  Anneal  Experiments 

It  is  instructive  at  this  point  to  consider  the  types  and  amounts  of 
radiation  damage  which  result  from  a typical  neutron  transmutation  irradiation. 
We  will  consider  the  following  sources  of  damage: 

• Fast  neutron  damage  for  a graphite  moderated  neutron  energy  spectrum 
(this  case  will  be  an  approximation  to  our  sample  irradiation  position). 

• Fast  neutron  damage  from  a fission  neutron  energy  spectrum  (this  case 
will  approximate  in  pile  irradiations  found  in  the  literature). 

• Gamma  damage  from  fission  gammas  from  the  reactor  pile. 

• Gamma  recoil  damage  (this  damage  occurs  as  a result  of  conservation  of 
momentum  when  a gamma  is  emitted  from  an  excited  Si  nucleus  after 
thermal  neutron  capture). 

• Beta  recoil  damage  (this  damage  occurs  as  a result  of  conservation  of 

momentum  when  the  beta  is  emitted  in  the  deexcitation  of  ^Si  ■+■  + 8~. 

We  will  neglect  damage  from  (n,  p),  (n,  a),  ana  ( y,  n)  etc.  , reactions 

because  of  their  low  cross  sections  for  neutron  energies  of  interest. 

Although  (n,  y)  reactions  with  boron  have  large  cross  sections,  the  boron 
concentration  in  our  samples  is  negligible  compared  to  the  silicon  concen- 
tration. Also,  we  estimate  that  the  damage  from  the  8 emitted  during 
31 

P production  will  produce  about  the  same  number  of  displacements  as  the 
recoil. 


for  Ed  ^ 2E  where  E is  the  average  kinetic  energy  of  the  recoiling 
R d R 

lattice  atom  and  E,  is  the  displacement  threshold  in  silicon  which  we 
d 22 

take  as  about  12  eV.  The  maximum  recoil  kinetic  energy  of  a target  atom 
struct  by  a neutron  of  energy  E^  is  given  by 


Ep(max) 


y 4e 


(23) 


when  A is  the  target  material  atomic  weight. 
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For  a fission  neutron  spectrum,  the  average  recoil  kinetic  energy  is 

v 23 
given  by 

- E (max) 

R (fission  neutron)  = 1 (24) 


while  for  a graphite  moderated  spectrum 
E 


(moderated  neutron)  = ER(max) 


In 

Ep(max)- 

L d J 

(25) 


We  shall  assume  that  the  average  neutron  energy  to  be  used  in  equation  (23) 
is  1.5  MeV.  Then  we  find  from  (23)»  (24),  and  (25)  that 

E (max)  = 4(1.5  MeV)  = 214  keV, 


R ~2T 

(fission  neutron)  = 214  keV 


107  keV 


and 


ER  (moderated  neutron)  = 1/2(214  keV) 


In 


214  x 10  eV 


= 10.9  keV. 


12  eV 

-J 

We  find,  using  (22),  that  the  number  of  displacements  per  neutron  is 
much  higher  for  a fission  spectrum  than  for  a moderated  spectrum,  i.e., 

107  keV 

V (fission)  = ' ~2~(~i2' ~ev ) = **^58  displacements /neutron  (26) 

and 

v (moderated)  = g^'^ev)^  = ^55  displacements /neutron.  (27) 

Furthermore,  in  our  bulk  pool  facility,  assuming  a Cd  ratio  of  about  30:1, 
the  fast  flux  is  about  (5  x 101:Ln/cm2/sec)/30  = 1.67  x 1010n/cm2/sec.  For 
an  in-core  irradiation,  the  Cd  ratio  is  about  1:1  and  a typical  fast  flux 
is,  therefore,  5 x 10  n/cm  /sec.  Using  these  fluxes  and  assuming  that  the 

displacement  cross  section  is  approximately  equal  to  the  scattering  cross 

23  -24  2 

section  , which  we  will  assume  is  3 x 10  cm  , we  find  using  (21)  and  (26) 
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(fast  neutron  in  core)  = (5  x 102t')(3  x 10  *^)(UU58)(5  x lO1^) 
dt 

= 3.35  x 10  displacements/cin  /sec.  (28^ 
Using  (21 ) and  (27  ) we  find  also  that 


dN. 

dt 


D pp  pL  10 

fast  neutron  in  pool)  = (5  x 10  )(3  x 10  )( 1+55)  (1.67  x 10'  ) 

12  3 
= l.lU  x 10  displacements/cin  /sec. 


(29) 


The  gamma  flux  radiation  damage  is  also  quite  different  since  in  core 

the  photon  flux  is  about  3 x lO^R/sec  while  in  our  bulk  pool  facility  it  is 
2 

about  3 x 10  R/sec.  The  cross  section,  o(y)»  for  displacements  is  dominated 

by  Compton  electron  displacements  and  is  nearly  independent  of  photon 

23  2l+ 

energy  in  silicon  for  typical  fission  gamma  energies.  We  will  assume 

-25  2 

that  this  cross  section  is  about  1 x 10  cm  for  100  keV  < E <2  MeV  and 

9 2 25  ^ 

that  1 R/sec  = 2 x 10  1 MeV  photons/cm  /sec.  Then  the  number  of 

displacements  is 


dN 

dt~ 


( gamma  in  pile)  = (5  x 10^)(l  x 10  2^)(3  x 10^) (2  x 10^)  (30) 

13  3 

= 3 x 10  displacements/cm  /sec 


and 


dN 

dt 


D (gamma  bulk  pool)  = (5  x 10^) (l  x 10  ^'(3  x 1C  )(2  x 10^) 


9 3 

= 3 x 10  displacements/cm  /sec. 


(31) 


The  various  silicon  isotopes  after  thermal  neutron  absorption  relax 
from  excited  states  by  prompt  gamma  emission.  This  imparts  a recoil  to  the 
silicon  isotope  in  order  to  conserve  momentum.  Since 


Ey 

= c E MV* 


Ed  = 1/2MV2 


1/2 


EY 

? 

MC 


An  average  over  all  silicon  isotopes  and  cross  sections  yields  an  average 

Ok  27 

recoil  kinetic  energy  of  780  eV.  ’ Therefore,  using  (22),  we  find 


— _ .780  eV 
V 2(l2  eV) 


32.5  displacements/thermal  neutron. 


I' 
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The  average  cross  section  for  these  processes  is  1.3  x 10  cm  and  the 
flux  of  interest  is  the  thermal  neutron  flux,  therefore 

( Y-recoil  in  core)  = (5  x 1022)(l.3  x 10-25)(5  x 1013)(32.5) 

dt  ^ 3 3 

= 1.06  x 10  displacements/cm  /sec  (32) 

and 

dN 

( y-recoil  in  pool)  = (5  x 1022)(1.3  x 10_25)(5  x 101;L)(32.5) 
clx> 

11  3 

= 1.06  x 10  displacements /cm  /sec.  (33) 


A similar  procedure  can  be  applied  to  beta  decay  recoil.  Assuming  the 

beta  energy  to  be  E„  = 1.5  MeV,  and  using  conservation  of  momentum,  we 
p 

find 

P.  /e62-  (m0C2)2 

c 

EL  = 1/2MV2  = 1/2 

R 

Therefore,  using  (22  ), 

33  2 eV 

V = 2(12  eV)  = 2.76  displacements/neutron  absorbed. 

Since  this  occurs  for  only  3%  of  the  total  silicon  atoms,  with  a cross 

-25  2 

section  of  1.1  x 10  cm  , 


= MV 


E 2 - (m  C2)2 

B o = 33.2  eV. 


MC 


and 


(6 


recoil  in  core)  = ( 0. 03) ( 5 x 1022)(l.l  x 10_25)(5  x 1013)(2.76) 
10  3 

= 2.3  x 10  displacements/cm  /sec  (3^) 


dt 


(B  recoil  in  pool)  = (0.03)(5  x 1022)(l.l  x 10-25)(5  x 1011 

8 3 

= 2.3  x 10  displacements/cm'"/sec. 


)(2.?6) 

(35) 


These  numbers  must  be  compared  with  the  rate  of  producing  phosphorus  atoms 
which  is 


dNp 

(in  core).  = 


(0.03)(5  x 1022)(1.1  x 10~25)(5  x 1013) 
8.25  x 109  [p]/cm3/ sec. 


(36) 
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and 


P ??  Pr»  11 

— • (in  pool)  = ( 0.03) ( 5 x 10  )(l.l  x 10  J){5  x 10  ) 

= 8.25  x 10T  [p]/cm3/sec.  (37) 

We  can  now  compare  the  number  of  displacements  per  phosphorus  atom 
produced  by  each  of  these  mechanisms.  This  information  is  summarized  in 
Table  11. 

Although  these  calculations  are  quite  speculative,  several  differences 
between  in-core  (low  Cd  ratio)  and  in-pool  (high  Cd  ratio)  irradiations 
become  apparent.  It  is  clear  that  the  fission  gamma  and  beta  recoil 
damage  is  point-defect  like.  The  annealing  characteristics  of  this  damage 
should  be  similar  to  1 Mev  electron  irradiation  data.  It  is  also  clear  that 
the  gamma  recoil  damage  is  a significant  fraction  of  the  radiation  damage 
for  in-pool  irradiations  but  is  relatively  less  significant  for  in-core 
irradiations . 

Chukichev  and  Vavilov  have  estimated  experimentally  that  the  number  of 

gamma  recoil  displacements  is  about  equal  to  the  number  of  fast  neutron 

27 

displacements  in  a heavy  water  moderated  reactor.  Our  ratio  of  fast 
neutron  to  gamma  recoil  displacements  is  of  the  order  of  10:1  for  an 
in-pool  position  with  a Cd  ratio  of  30:1.  Since  the  Cd  ratio  in  a heavy 
water  reactor  is  between  100:1  to  1000:1,  our  calculations  are  in  reasonable 
agreement  with  the  experimental  data  in  Ref.  27. 

It  is  not  clear  whether  to  treat  the  gamma  recoil  radiation  damage  as 
point  defect  like  or  cluster  like.  The  average  number  of  displacements  per 
neutron  is  rather  large  (32.5  displacements /thermal  neutrcr ' and 

suggests  at  least  the  production  of  small  clusters.  The  actual  number  of 
displacements  per  thermal  neutron  absorbed,  however,  range  from  about  10 
or  20  up  to  80  to  160.  There  is,  therefore,  a continuous  distribution  of 
cluster  sizes  which  approach  the  point  defect  distributions  at  the  lower 
end.  The  same  cluster  size  arguments  can  be  made  about  fast  neutron  damage 
from  a highly  moderated  energy  spectrum. 

We  have  considered  gamma  recoil  damage  both  as  cluster  like  and  point 
defect  like  in  Table  11,  to  calculate  the  ratio  of  cluster  like  displacements 
to  point  defect  like  displacements  for  our  two  cases  of  interest.  It  seems 
likely  that  the  in-core  irradiations  deposit  300  to  3000  more 
displacements  in  clusters  than  in  point  defects  as  compared  to  in-pool 
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TABLE  11.  Number  of  Displacements  per  Phosphorus  Produced 


Fast  Neutron 

U.06  x 10^ 

L 

1.38  x 10 

Fission  Gamma 

1.0 

36.  It 

Gamma  Recoil 

1.29  x 103 

1.29  x 103 

Beta  Recoil 

2.76 

2.76 

Total  Disp./  [ P] 

4.06  x 10^ 

I4 

1.51  x 10 

Cluster/Pt.  Defect  Ratio: 

(if  gamma  recoil  = cluster) 

31^0 

10.  u 

(If  gamma  recoil  = pt.  defect) 

1.08  x 106 

385 

Ratio  of  Core  to  Pool  Cluster  to  Pt.  Defect  Ratios: 
(Gamma  recoil  = cluster)  302 
(Gamma  recoil  = pt.  defect)  2805 


I 

I 

I I 


I 


i 

J 

i 

i 

i 

I 

f 

[ 


I 
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irradiations.  Such  a large  difference  in  the  defect  spectrum  is  very  likely 
to  produce  differences  in  annealing  characteristics  for  the  two  cases. 

We  will  discuss  some  of  these  differences  in  the  following  discussion  on 
annealing. 

Table  12  is  a summary  of  ali  the  irradiations  performed  thus  far, 

the  sample  number,  the  fluence,  the  phosphorus  concentration 

calculated  from  the  cross  section,  and  initial  resistivity  assuming  a 

typical  mobility,  and  the  purpose  of  the  irradiation.  In  this  table, 

•fp  is  the  fluence  required  to  take  the  particular  p-type  sample  from  its 

initial  resistivity  to  intrinsic  (or  exact  compensation  of  the  residual 

boron  by  the  phosphorus  added) . An  inspection  of  this  table  shows  immediately 

the  failure  to  control  the  resistivity  for  partial  compensation  on  the  p- 

type  side  of  intrinsic  (i.e.  <J>  < <t^,)  due  to  premature  type  conversion. 

Figure  29  shows  this  premature  type  conversion  effect.  This  figure 

compares  an  isochronal  anneal  in  argon  of  a wafer  of  Rockwell  float  zone, 

MURR  1,  with  the  annealing  in  air  of  a piece  of  Monsanto  Czochralski,  and 

2 8 2Q 

vacuum  anneals  of  float  zone  and  Czochralski  by  Kharchenko  et  al.  ’ y 

The  solid  bar  is  the  expected  p-type  resistivity  for  the  Rockwell  float 

zone  after  the  final  anneal  on  the  basis  of  the  number  of  donors  added 

[0 . 0U2  ppb].  Instead  the  final  resistivity  is  about  HOO  Q-cm  n-type  (the 

type  conversion  probably  occurred  between  600  and  800°C)  a gain  of  extra 

13  3 

donors  of  0.282  ppb  (l.Ul  x 10  atoms/cnr)  beyond  that  added  by  neutron 
doping  [0.0U2  ppb]. 

Also  shown  on  this  figure  is  the  isochronal  annealing  of  a second 

unirradiated  Rockwell  sample  from  the  same  wafer.  It  should  be  noted 

that  the  unirradiated  sample  shows  no  adverse  effects  of  the  annealing 

and  that  the  resistivity  of  the  irradiated  sample  MURR  1 drops  below  that 

of  the  unirradiated  sample  at  300°C. 

The  heavily  irradiated  float  zone  sample  by  Kharchenko  shows  no 

type  conversion  peaks,  i.e.,  no  rise  in  resistivity  toward  intrinsic 

(230,000  fi-cm).  This  is  confirmed  by  their  Hall  measurements  directly 

29 

after  irradiation  but  before  annealing.  It  should  be  noted  that  a 

drafting  error  has  apparently  been  made  in  Figure  1 of  Ref.  29  since  the 

13  -3 

final  values  of  n = 5 x 10  cm  and  p = 1 fl-cm  imply  an  electron 

2 

mobility  of  125,000  cm  /V-sec,  some  two  orders  of  magnitude  too  high! 
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TABLE  12.  Summary  of  Irradiations 


Resistivity 
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29.  A comparison  of  isochronal  annealing  of 
transmutation  doped  float  zone  and 
Czochralski  silicon. 


19  2 

Their  sample  was  irradiated  to  a fluence  of  about  2.68  x 10  n/cm 

l6  2 

compared  to  the  Rockwell  sample  fluence  of  1.30  x 10  n/cm  . Therefore, 
their  sample  experienced  a fluence  of  some  2000  times  that  of  most  of  the 
samples  listed  in  Table  12.  Furthermore,  their  sample  was  irradiated  in 
core  while  ours  was  irradiated  in  pool  (see  Table  11).  Because  of  the 
fluence  differences  and  position  differences,  it  is  estimated  that  the 
number  of  displacements  going  into  clusters  is  6 x lO'’  to  6 x 10^  times 
greater  than  for  our  samples.  The  number  of  point  defects  is  about  2000 
times  greater. 

It  is  instructive  to  compare  the  radiation  damage  in  the  float  zone 

sample  by  Kharchenko,  et  al.  with  the  damage  produced  by  a P ion  implanted 

sample.  The  neutron  fluence  for  a 1 fl-cm  n-type  final  resistivity  and 

1000  fi-cm  p-type  starting  resistivity  can  be  calculated  from  (7)  to  be 
19  2 

2.68  x 10  n/cm  . The  number  of  displacements  is  then,  using  $ = <J>t, 
and  equation  (21), 

ND  = (5  x 1022)(3  x 10~2l+) (2.68  x 1019)(l*l+58) 

22  , 3 

= 1.79  x 10  displacements /cm  . (37) 

22  3 

Since  silicon  contains  only  5 x 10  atoms/cm  , it  is  clear  that  the 
radiation  damage  in  the  Kharchenko  sample  is  approaching  the  amorphorus 
condition. 

j To  produce  an  amorphorus  layer  by  P+  ion  implantation,  a dose  of  about 

5 x 10  ions/cm  at  Uo  keV  is  required.  The  number  of  displacements  per 

, ion  can  be  calculated  from  (22)  to  yield 

_ E _ 1*0  keV  . ,,,  . _3  .. 

v ~ 2E^  " 2Tl2  eV)  = 1,666  x 10  cusp.  /ion. 

Since  the  depth  of  a Uo  keV  P+  implant  is  about  2 x 10_?  cm,  because  of  en- 
hanced diffusion  (see  Figure  2.27,  p.  56  in  Ref.  30),  the  number  of 
displacements  is 

f 

lU  2 3 

= ( 5 x 10  ions /cm  )( 1.667  x 10  displacements /ion) 

2 x 10  9 cm 

22  2 

= 4 x 10  displacements/cm  . (38) 

We  can  conclude  from  a comparison  of  (37  ) and  ( 38  ) that  the  Kharchenko 
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float  zone  sample  is  nearly  amorphorus  as  a result  of  fast  neutron  damage. 

This  condition  has  interesting  consequences.  In  ion-implanted  silicon, 

annealing  up  to  temperatures  of  about  800°C  are  required  for  complete 

electrical  activity  of  the  implanted  ion  in  cases  where  an  amorphorus  layer 
31 

is  not  formed.  The  recovery  of  the  electrical  activity  of  an  amorphorus 
layer  is  completed  by  about  600°C  and  is  closely  related  to  the  formation 

of  crystalline  films  of  Si  during  vacuum  evaporation  and  to  epitaxial 

..  33 
regrowth. 

Fig  ire  29  also  shows,  for  comparison,  the  isochronal  annealing  of  two 

10  3 

Czochralski  silicon  samples  ([o]  'u  10  cm  J).  The  data  for  the  Kharchenko 
sample  is  taken  from  Ref.  29.  The  reverse  annealing  peaks  are  evident  in 
both  Czochralski  samples  and  are  apparently  related  to  the  high  oxygen 
concentration  in  these  samples  and  not  to  neutron  fluence. 

Evidence  for  this  point  of  view  is  shown  in  Figure  30.  These  float 
zone  samples  were  irradiated  to  a fluence  sufficient  to  produce  'u  100  ft-cm 
n-type.  The  expected  final  resistivity  values  calculated  from  the  starting 
resistivity  are  shown  on  the  figure.  The  Topsil  sample  came  to  within  1? 
of  the  expected  donor  concentration  while  the  Rockwell  sample  was  only 
slightly  worse.  It  should  be  noted  that  the  fluences  for  these  two  samples 
are  similar  to  the  fluences  for  the  Czochralski  samples  in  the  previous 
figure,  however,  the  annealing  peaks  and  dips  between  200  and  500°C  do  not 
appear  in  the  float  zone  material  at  this  fluence.  Furthermore,  dislocation 
density  does  not  appear  to  have  a direct  correlation  with  the  lack  of  the 
reverse  annealing  peaks  shown  in  Figure  29  since  the  Rockwell  sample  in 
Figure  30  is  highly  dislocated  while  the  Topsil  sample  is  not.  We  conclude 
then,  on  the  basis  of  circumstantial  evidence,  that  the  reverse  annealing 
peaks  are  a result  of  the  concentration  of  oxygen  relative  to  the  level  of 

TO 

phorphorus.  (The  oxygen  concentration  is  thought  to  be  between  10  to 
15  -3 

10  cm  in  float  zone  silicon). 

Kharchenko  et  al.  also  state  that  their  Czochralski  samples  were  p- 

29 

type  up  to  600°C  and  n-type  thereafter.  (Their  Czochralski  samples  are 
apparently  irradiated  in-pool).  We  conclude  from  our  data  in  Figure  30 
that  electrical  activity  of  the  phosphorus  begins  at  about  600°C  in 
agreement  with  their  conclusions.  The  large  reverse  annealing  peak  from 


RESISTIVITY  (fl-CM) 


Figure  30.  Isochronal  annealing  of  two  float  zone  samples 
irradiated  to  produce  ^100  .1-cm  N-type. 
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600-800°C  is  primarily  associated  with  the  electrical  activation  of 
phosphorus.  The  Rockwell  sample  resistivity  is  intrinsic  at  700°C  while 
the  intrinsic  peak  in  the  Topsil  sample  has  apparently  been  missed.  It 
can  further  be  concluded  from  the  data  in  Figure  30  that  nearly  all  the 
electrically  active  radiation  damage  has  annealed  at  temperatures  belcw 
600°C.  This  contention  can  be  justified  as  follows. 

We  can  calculate  the  number  of  electrically  active  donors  as  a function 
of  isochronal  annealing  temperature  assuming  that  the  mobility  has  completely 
recovered  before  600°C  and  that  the  reverse  annealing  peak  between  600 
and  800° C is  the  result  of  a type  conversion  from  p-type  to  n-type  due 
to  the  electrical  activation  of  phosphorus.  (The  thermal  probe  type  as 
a function  of  annealing  temperature  is  shown  on  the  figure).  On  the  basis 
of  this  assumption  we  can  calculate  the  fraction  of  phosphorus  which  is  not 
electrically  active  (N/No)  as  a function  of  annealing  temperature.  We 
assume  that  this  fraction  is  unity  at  the  resistivity  minimum.  Then  if  the 
electrical  activation  is  governed  by  first  order  annealing  kinetics. 


d(N/N  ) 

o_ 

at 


-k(t)(n/n  ) 

O 


— E/kT 

where  K(T)  = N e is  the  temperature  dependent  rate  at  which  the 

fraction  of  phosphorus  not  electrically  active  decreases  at  the  various 
annealing  temperatures.  From  the  above 


in  (N  /N)  = (N  t)e'E/k,r 

o o 

where  t is  the  annealing  time  at  temperature  T.  For  equal  annealing  times 


at  various  temperatures,  a plot  of  In  N^/N  vs.  1/T  yields  the  activation 
energy  for  the  process  while  the  intercept  yields  the  defect  jump  frequency. 


The  result  of  this  first  order  analysis  is  shown  in  Figure  31.  It 
is  clear  that  the  first  order  annealing  assumption  is  reasonable  and  that 
the  first  order  annealing  occurs  throughout  the  type  conversion  peak  in 
Figure  30.  To  have  obtained  this  result,  the  two  assumptions  mentioned 
previously  must  be  valid,  i.e.,  e lack  of  electrically  active  damage  at 
temperatures  above  600°  and  the  complete  recovery  of  carrier  mobility 
before  600°C. 

The  first  order  annealing  analysis  suggests  that  the  activation  energy 
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Figure  31.  Semilog  plot  of  the  log  of  the  reciprocal 

of  the  fraction  of  phosphorus  not  electrically 
activ'  vs.  reciprocal  annealing  temperature 
for  the  Topsil  sample  shown  in  Figure  30. 
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for  the  phosphorus  to  become  electrically  active  is  O.96  + 0.05  eV  while 


2 -1 

the  jump  frequency  is  about  2 x 10  sec 


This  activation  energy  is 


remarkably  close  to  the  values  given  in  the  literature  for  the  reorientation 

of  the  vacancy  phosphorus  E-center  (E  = 0.93  +_  0.05  eV,  = 1.6  x 10  sec  “)J 

and  for  the  electrical  disappearance  of  a deep  acceptor  (E^-O.MeV  thought 

to  be  the  E-center  (E  = Q.91*  eV,  V = lO^sec-1).^ 

Watkins  and  Corbett  suggest  that  the  lower  Jump  frequency  in  Ref.  35 

indicates  the  diffusion  of  the  E-center  through  the  lattice  over  great 

distances  (10~^/10^  = 107  jumps  before  dissociation  at  a vacancy  sink). 

They  point  out  that  an  E-center  reorientation  followed  by  a vacancy  phosphorus 

interchange  is  equivalent  to  an  E-center  jump  through  the  lattice.  If  we 

13  2 1C 

make  this  same  assumption,  then  we  are  observing  1C)  /2  x 10  = 5 x 10 

jumps  before  annihilation!  For  a random  walk  of  5 x 10^  jumps  of  a distance 

O 

of  a,  2A,  the  mean  distance  traveled  would  be 


5 x 10 


(2  x 10  ^cm)  = U. 5 x 10  Jcm. 


This  is  of  the  order  of  the  smallest  sample  dimension  and  would  suggest 
vacancy  anihilation  at  the  sample  surface.  But  this  can  not  be  the  correct 
mechanism  for  E-center  annihilation  because  it  would  deposit  all  the 
phosphorus  at  the  sample  surface.  We  know,  however,  that  phosphorus  is 
uniformly  distributed  in  transmutation  doped  silicon.  We  must  conclude, 
then,  that  free  migration  of  E-centers  in  NTD-Si  does  not  occur  for  some 
reason.  One  possibility  is  that  the  E-center  becomes  pinned  in  some  way 
to  its  lattice  position  by  the  clustering  with  other  defects.  Since  there 
are  approximately  three  displacements  per  beta  emission  when  the  phosphorus 
atom  is  created  by  transmutation,  it  is  possible  that  a di-interstitial  or 
di-vacancy.  E-center  complex  is  formed  which  remains  stable  up  to  the 
higher  annealing  temperatures  observed  here. 

It  should  be  noted  in  Figure  30  that  at  the  minimum  resistivity  near 

o 1^3 

600°C  in  the  Topsil  sample,  the  hole  concentration  is  about  6 x 10^  cm~J. 

Since  the  initial  hole  concentration  was  about  2 or  3 x 10  cm  \ we  can 
safely  infer  that  a very  large  acceptor  concentration  has  been  created 
as  a result  of  irradiation  and  annealing  to  550°C  (an  alternate  interpretation 
would  be  the  destruction  of  6 x 10x  cm  donors  but  this  would  imply  a 
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compensation  ratio  of  0.9996  before  irradiation  which  is  very  unlikely). 

It  will  be  shown  shortly  that  similar  acceptor  defect  concentrations 
have  been  observed  in  all  float  zone  studied  to  date  and  that  this  defect 
concentration  is  independent  of  irradiation  fluence  over  a fluence  range 
of  50:1.  We  conclude,  therefore,  that  the  magnitude  of  this  concentration 
is  dependent  on  the  properties  of  the  silicon  before  irradiation.  The 

most  likely  impurity  to  be  found  in  highly  zone  refined  silicon  in  a 

1 3 1 ^4  3 

concentration  range  of  3 x 10  cm  (Rockwell,  Wacker)  to  5 x 10x  cm- 

(Topsil  and  T.I.)  is  oxygen.  We,  therefore,  suggest  that  these  acceptor 
defects  are  single  or  multiple  defect-oxygen  complexes,  i.e.,  the  oxygen 
impurity  concentration  is  decorated  as  a result  of  the  process  of  radiation 
damage  and  annealing  to  550° C. 

Figure  32  shows  the  recovery  of  resistivity  and  minority  carrier 
lifetime  on  a Rockwell  sample  l"  dia  x 1"  long.  The  fluence  for  this 
irradiation  was  50%  of  that  necessary  to  produce  type  conversion.  The 
annealing  shows  several  unusual  features  which  are  found  only  in  the 
Rockwell  samples  and  which  could  possibly  be  related  to  the  dislocation 
density.  The  resistivity  between  200°C  ans  800°C  is  multivalued.  Shortly 
after  removal  from  the  annealing  furnace  and  quenching  to  room  temperature 
in  trichloroethylene,  the  resistivity  is  a minimum  but  rises  to  the  second 
higher  value  while  storing  the  sample  in  the  dark.  Upon  illumination  from 
a tungsten  light,  the  resistivity  returns  to  its  t = o value.  The 
resistivity  can  be  cycled  any  number  of  times.  The  reciprocal  of  the 
minority  carrier  lifetime  follows  the  resistivity  in  general.  Very  heavy 
trapping  is  observed  in  the  same  region  as  the  largest  change  in  resis- 
tivity. This  sample  type-converts  at  about  800°C.  It  should  al.  o be 
noted  that  a plateau  in  the  recovery  around  300°C  is  similar  to  the 
Czochralski  sample  annealing  seen  previously. 

In  an  attempt  to  determine  if  sample  contamination  is  responsible 
for  the  unexpected  type  conversion,  a series  of  isochronal  anneals  were 

_7 

performed  in  argon  and  a vacuum  of  5 x 10  torr.  The  fluences  for  these 

irradiations  were  in  the  range  of  $ = 0.36  <J>  to  $ = 0.56  $ where  $ is 

c c c 

the  fluence  required  for  exact  compensation.  The  results  are  shown  in 
Figures  33  through  39-  A comparison  of  a pair  of  vacuum  and  argon  anneals 
are  shown  in  Figures  Uo  and  Ul.  It  can  be  seen  that  all  of  these 
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Figure  32.  Isochronal  annealing  of  resistivity  and 

minority  carrier  lifetime  of  Rockwell  float 
zone  (argon  anneal)  - bulk  sample. 
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Figure  33.  Argon  isochronal  anneal  of  resistivity  of 
Rockwell  float  zone  wafer. 
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Figure  35.  Argon  isochronal  anneal  of  T.  I.  Lopex  wafer. 
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Figure  37.  Vacuum  isochronal  anneal  of  T.  I.  - Lopex  wafer. 
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Figure  39.  Vacuum  isochronal  anneal  of  Topsil  wafer. 
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samples  type  convert  at  temperatures  between  600  and  900°C.  The  expe^tei 
resistivities  calculated  from  the  concentration  of  phosphorus  added  is 
also  shown  on  this  series  of  figures. 

Note  that  the  minimum  occurs  at  the  same  temperature  in  all  of  these 
anneals  and  that  it  is  identical  to  the  heavily  irradiated  samples  shown 
in  Figure  30.  This  minimum  resistivity  varies  between  20C  to  2000  h-cm 
in  all  these  samples  Justify  our  estimates  of  the  concentration  of 
acceptor  defects  discussed  previously. 

An  unusual  feature  seen  in  most  of  these  anneals  cure  the  Czochralskl 
like  reverse  anneals  at  temperatures  below  600°C.  These  reverse  anneals 
are  most  pronounced  in  those  samples  which  show  the  lowest  minimum^  at 
550  to  600°C.  This  is  evident  in  Figures  3**  and  35.  They  are  much  less 
pronounced  in  the  Rockwell  and  Wacker  silicon  which  we  suspect  indicates 
about  an  order  of  magnitude  less  oxygen  in  these  samples . 

The  .inal  resistivities  are  completely  uncontrollable  for  these  light 
irradiations.  Although  it  is  tempting  to  suspect  sample  contamination 
during  the  anneal,  we  have  been  unable  to  detect  any  contamination  by 
neutron  activation  analysis.  Furthermore,  Figure  Ul  suggests  that  there 
is  very  little  difference  between  argon  and  vacuum  annealing.  Figure  31* 
shows  a type  conversion  at  600°C.  It  is  difficult  to  imagine  diffusing 
3 x 10  atoms/cm  into  this  sample  from  surface  contamination  at  600°C. 

We  have  seen  evidence  of  surface  contamination  for  some  anneals  above 
950°C,  however,  impurity  diffusion  is  much  more  likely  at  this  temperature. 

Figure  U2  shows  isochronal  annealing  in  argon  for  an  irradiation 
of  <!>  = 1.02  $c,  i.e.  2%  overcompensation.  The  low  resistivities  at 
temperatures  below  300°C  were  the  result  of  remounting  this  sample  with 
wax  on  a hot  plate  after  irradiation  for  relapping  after  etching.  It 
is  seen  that  the  final  resistivity  is  very  close  to  the  expected  value 
except  for  the  last  annealing  point  which  may  be  surface  contamination. 
Because  of  slight  overcompensation,  type  conversion  is  expected  in  this 
sample.  The  etching  and  lapping  may  have  removed  a sodium  surface 

v 

contamination,  however,  the  same  procedure  on  other  samples  has  made  no 
improvement  in  our  ability  to  prevent  type  conversion. 
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Perhaps  the  additional  phosphorus  concentration,  the  extra  radiation 
:ana**e  or  the  final  charge  type  has  suppressed  the  unknown  donor  vhi  eh  '■  -se: 
type  conversion  in  some  way. 

Tn  the  following  table,  we  summarize  some  of  the  prominent  annealing 
features  common  to  all  samples  investigated  tc  date.  The  parameters 
listed  below  are  tne  two  "Czochra  ski-like"  peaks,  the  first  ar.nealir.ff 
temperature  at  which  the  resistivity  equals  the  before  irradiation 
resistivity,  the  minimum  resistivity  temperature  and  the  ratio  of  the 
minimum  resistivity  to  the  before  irradiation  resistivity. 


V.  F?  OPESS  AND  ACCOMPLISHMENTS 

The  object  of  this  research  program  is  to  develop  the  technique  of 
neu'r.-:.  transmutation  doping  in  silicon  to  the  point  that  the  compensation 
rati  can  be  controlled  to  a very  high  precision.  Significant  success  in 
this  task  will  al-ov  much  better  control  over  the  boron  compensation  in  silicon 
lsed  for  extrinsic  I?  detection.  This  in  turn  should  provide  a better  control 
over  detector  detectivity  and  temperature  dependence. 

V..  ugh  considerable  unexpected  difficulty  r.as  been  encountered  ir. 

• -rolling  the  inpurity  concentrations  in  silicon  with  any  degree  of 
■tair.ty  fceiov  0.1  ppb,  we  feel  that  real  progress  has  beer,  made  in  achieving 

• gca  . It  has  also  become  evident  that  a better  understanding  of  the 
nacla*  . •:  damage , the  oxygen  concentration  of  the  samples  and  the  interaction 
betweer  the  defects  and  residual  inpurities  in  silicon  before  irradiation  are 
nee  5e  . 1 “ier  to  establish  control  ever  the  entire  doping  process  ir. 
concentrations  below  0.1  ppb. 

We  believe  that  neutron  transmutation  doping  has  removed  one  of  the  major 
-ifficulties  ir.  achieving  exact  compensation  in  silicon,  i.e.,  the  precise 
control  over  phosphorus  concentration  which  is  improved  by  a factor  of  at  least 
50:1  over  conventional  doping  techniques.  Unfortunately,  other  residual 
impurity  problems  have  bee.,  detected  for  the  first  time  as  a result  cf  achieving 
this  control.  Because  of  these  difficulties,  we  are  behind  schedule  ir. 
compensating  the  boron  in  the  presence  cf  a very  high  Ga  or  In  concentration. 

We  fee_,  . -.ever,  that  a much  better  -understanding  of  the  limitations  involved 
in  this  problem  have  resulted  from  our  concerted  effor*  or.  undoped  material 
where  impurity  concentrations  car.  be  determined  with  a much  higher  precision. 

We  also  believe  tr.at  we  now  neve  evidence  tc  support  the  contention  that 
compensation  control  ir.  r.-type  material  is  possible. 

The  following  is  ar.  abbreviated  list  of  the  significant  accomplishments 
achieved  during  this  program. 

• Construction  of  the  Bulk  Po^l  Variable  Flux  Facility  and  demonstration 
of  fluer.ee  reproducibility  to  +_  0.^. 

• Demonstration  of  dopir.c  concentration  control  to  ♦_  12  for  ICO  ft-cm 
material. 

• So  significant  impur'ty  problems  detectable  by  neutron  activation 
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analysis  a*.wr  tc  concentrations  of  tetter  than  0.1  ppb  with  the  exception 

13  3 

of  sodium  ir.  concentrations  of  1 to  6 x 10  atoms/cm  . 

• Pre- irradiation  materials  analysis  of  p,  p,  and  T and  activation 
analysis  indicate  that  silicon  from  all  sources  are  roughly  compariable  in 
quality.  The  largest  exception  is  in  minority  carrier  lifetime,  t,  which 
is  not  as  large  as  advertised. 

• Theoretical  analysis  of  the  behavior  on  n,  r.,j,  y,  and  p versus  f-uence 
has  been  developed  to  evaluate  results  of  neutron  doping  compensation  of 
undoped  materials. 

• The  theoretical  errs  section  for  phosphorus  production  has  been 
verified  to  about  +_  1^. 

• Estimates  have  been  made  on  the  number  of  displacements  as  a function 
of  fast  to  thermal  neutron  ratio.  Importance  of  gamma  recoil  damage 
calculations  have  been  verified. 

• Several  differences  have  been  discovered  between  our  annealing  results 
and  those  found  in  the  literature.  These  have  been  explained  in  terms  of 
fast  to  thermal  neutron  ratio  and  carrier  type  immediately  after  irradiation 
but  before  annealing. 

• Have  related  in  core  1 fJ-cm  radiation  and  subsequent  annealing  results 
discussed  in  literature  to  a radiation  damage  produced  amorphous  condition 
which  alters  annealing  kinetics  drastically. 

• The  activation  energy  and  frequency  factor  for  the  electrical  activation 
of  transmutation  produced  phosphorus  during  annealing  has  been  determined 
for  the  first  time.  This  analysis  suggests  that  vacancy-phosphorus  pairs 
(E-centers)  are  related  to  the  process  of  recovery  of  electrical  activity. 

• "Czochralski-like"  annealing  has  been  observed  in  lightly  irradiated 
silicon  ($  t 0.5  $r)  which  suggests  that  clustering  of  radiation  defects 
with  residual  oxygen  is  a problem  which  must  be  solved  in  float  zone 
doped  to  concentrations  below  0.1  ppb. 

13 

• The  production  of  an  acceptor  defect  in  concentrations  between  5 x 10 

lU  -3 

to  5 x 10  cm  by  irradiation  and  annealing  has  been  demonstrated  to  be 
independent  of  fluence.  This  strongly  suggests  that  this  acceptor  is 
related  to  an  intrinsic  residual  impurity  concentration  in  the  silicon  before 
irradiation.  We  propose  that  oxygen,  which  has  been  decorated  with  defects 
to  make  it  electrically  active,  is  the  source  of  these  acceptors.  This 
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result  is  exciting  since  it  suggests  that  oxygen  detection  in  silicon 
can  easily  be  extended  to  lower  than  present  day  detection  limits  by 
several  orders  of  magnitude. 

We  have  been  able  to  produce  a slightly  overcompensated  sample  in 
ax  least  one  case  where  the  ratio  of  donor  to  acceptor  concentration  is 
i.e.,  a control  of  compensation  ratio  close  to  2%  of  exac-  compensator. . 
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